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Abstract. We generalize the Newtonian n-body problem to spaces of cur- 
vature K = constant, and study the motion in the 2-dimcnsional case. For 
K > 0, the equations of motion encounter non-collision singularities, which 
occur when two bodies are antipodal. This phenomenon leads, on one hand, 
to hybrid solution singularities for as few as 3 bodies, whose corresponding 
orbits end up in a collision-antipodal configuration in finite time; on the 
other hand, it produces non-singularity collisions, characterized by finite ve- 
locities and forces at the collision instant. We also point out the existence 
of several classes of relative equilibria, including the hyperbolic rotations 
for K < 0. In the end, we prove Saari's conjecture when the bodies are on a 
geodesic that rotates elliptically or hyperbolically. We also emphasize that 
fixed points arc specific to the case k > 0, hyperbolic relative equilibria to 
K < 0, and Lagrangian orbits of arbitrary masses to k = — results that 
provide new criteria towards understanding the large-scale geometry of the 
physical space. 
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1. Introduction 

The goal of this paper is to extend the Newtonian n-body problem of celes- 
tial mechanics to spaces of constant curvature. Though attempts of this kind 
existed for two bodies in the 19th century, they faded away after the birth 
of special and general relativity, to be resurrected several decades later, but 
only in the case n = 2. As we will further argue, the topic we are opening 
here is important for understanding particle dynamics in spaces other than 
Euclidean, for shedding some new light on the classical case, and perhaps 
helping us understand the nature of the physical space. 

1.1. History of the problem. The first researcher who took the idea of 
gravitation beyond R'^ was Nikolai Lobachevsky. In 1835, he proposed a Kepler 
problem in the 3-dimensional hyperbolic space, H^, by defining an attractive 
force proportional to the inverse area of the 2-dimensional sphere of the same 
radius as the distance between bodies, [52]. Independently of him, and at 
about the same time, Janos Bolyai came up with a similar idea, [1]. 

These co- discoverers of the first non-Euclidean geometry had no followers in 
their pre-relativistic attempts until 1860, when Paul Joseph Serrelo extended 
the gravitational force to the sphere S^ and solved the corresponding Kepler 
problem, [55]. Ten years later, Ernst Schering revisited Lobachevsky's law 
for which he obtained an analytic expression given by the cotangent potential 
we study in this paper, [62]. Schering also wrote that Lejeune Dirichlet had 
told some friends to have dealt with the same problem during his last years 
in BerlinG, [63]. In 1873, Rudolph Lipschitz considered the same problem in 
S^, but defined a potential proportional to l/sin(r/i?), where r denotes the 
distance between bodies and R is the curvature radius, [SI]. He obtained the 
general solution of this problem in terms of elliptic functions, but his failure 
to provide an explicit formula invited new approaches. 

In 1885, Wilhelm Killing adapted Lobachevsky's idea to S^ and defined an 
extension of the Newtonian force given by the inverse area of a 2-dimensional 
sphere (in the spirit of Schering), for which he proved a generalization of 
Kepler's three laws, [H]. Another contributor was Heinrich Liebmannjfl. In 
1902, he showed that the orbits of the two-body problem are conies in S'^ 
and H^ and generalized Kepler's three laws to k 7^ 0, [1?]. One year later, 
Liebmann proved S^- and H^-analogues of Bertrand's theorem, [3], [TB], which 
states that there exist only two analytic central potentials in the Euclidean 



Paul Joseph Serret (1827-1898) should not be confused with another French mathe- 
matician, Joseph Alfred Serret (1819-1885), known for the Frenet-Serret formulas of vector 
calculus. 

This must have happened around 1852, as claimed by Rudolph Lipschitz, |50| . 
'^ Although he signed his works as Heinrich Liebmann, his full name was Karl Otto Hein- 
rich Liebmann (1874-1939). He did most of his work in Heidelberg and Munich. 
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space for which all bounded orbits are closed, [IE]. He also summed up his 
results in a book published in 1905, [lU] . 

Unfortunately, this direction of research was neglected in the decades fol- 
lowing the birth of special and general relativity. Starting with 1940, however, 
Erwin Schrodinger developed a quantum-mechanical analogue of the Kepler 
problem in S^, [6l]. Schrodinger used the same cotangent potential of Scher- 
ing and Liebmann, which he deemed to be the natural extension of Newton's 
law to the spherqj. Further results in this direction were obtained by Leopold 
Infeld, [3S], [H]- In 1945, Infeld and his student Alfred Schild extended this 
problem to spaces of constant negative curvature using a potential given by 
the hyperbolic cotangent of the distance. A list of the above-mentioned works 
also appears in [66], except for Serret's book, [65]. A bibliography of works 
about mechanical problems in spaces of constant curvature is given in [69] . 

Several members of the Russian school of celestial mechanics, including Va- 
leri Kozlov and Alexander Harin, [13], [15], Alexey Borisov, Ivan Mamaev, and 
Alexander Kilin, [5], [6], [7], [8], [39], Alexey Shchepetilov, [67], [68], [69], and 
Tatiana Vozmischeva, [H] , revisited the idea of the cotangent potential for the 
2-body problem and considered related problems in spaces of constant curva- 
ture starting with the 1990s. The main reason for which Kozlov and Harin 
supported this approach was mathematical. They pointed out, as Schering, 
Liebmann, Schrodinger, Infeld, and others had insisted earlier, that (i) the 
classical one-body problem satisfies Laplace's equation (i.e. the potential is a 
harmonic function), which also means that the equations of the problem are 
equivalent with those of the harmonic oscillator; (ii) its potential generates a 
central field in which all bounded orbits are closed — according to Bertrand's 
theorem. Then they showed that the cotangent potential is the only one 
that satisfies these properties in spaces of constant curvature and has, at the 
same time, meaning in celestial mechanics. The results they obtained support 
the idea that the cotangent potential is, so far, the best extension found for 
the Newtonian potential to spaces of nonzero constant curvature. Our paper 
brings new arguments that support this view. 

The latest contribution to the case n = 2 belongs to Jose Cariiiena, Manuel 
Ranada, and Mariano Santander, who provided a unified approach in the 
framework of differential geometry, emphasizing the dynamics of the cotangent 
potential in S^ and H^, [S] (see also [TD], [33] )• They also proved that, in this 
unified context, the conic orbits known in Euclidean space extend naturally 
to spaces of constant curvature, in agreement with the results obtained by 
Liebmann, 1661. 



^"The correct form of [the] potential (corresponding to 1/r of the flat space) is known to 
be cot X," [M], P- 14. 
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1.2. Relativistic n-body problems. Before trying to approach this prob- 
lem with contemporary tools, we were compelled to ask why the direction of 
research proposed by Lobachevsky was neglected after the birth of relativity. 
Perhaps this phenomenon occurred because relativity hoped not only to an- 
swer the questions this research direction had asked, but also to regard them 
from a better perspective than classical mechanics, whose days seemed to be 
numbered. Things, however, didn't turn out this way. Research on the classi- 
cal Newtonian n-body problem continued and even flourished in the decades 
to come, and the work on the case n = 2 in spaces of constant curvature was 
revived after several decades. But how did relativity fare with respect to this 
fundamental problem of any gravitational theory? 

Although the most important success of relativity was in cosmology and 
related fields, there were attempts to discretize Einstein's equations and define 
an n-body problem. Reamrkable in this direction were the contributions of 
Jean Chazy, [I3|, TuUio Levi-Civita, [S], gS], Arthur Eddington, [27], Albert 
Einstein, Leopold Infelco, and Banesh Hoffmann, [28], and Vladimir Fock, [30] . 
Subsequent efforts led to refined post-Newtonian approximations (see, e.g., 
[15] . [TB] , [17]). which prove useful in practice, from understanding the motion 
of artificial satellites — a field with applications in geodesy and geophysics — to 
using the Global Positioning System (GPS), [TB] . 

But the equations of the n-body problem derived from relativity prove com- 
plicated even for n = 2, and they are not prone to analytical studies similar to 
the ones done in the classical case. This is probably the reason why the need 
of some simpler equations revived the research on the motion of two bodies in 
spaces of constant curvature. 

Nobody, however, considered the general n-body problenu for n > 3. The 
lack of developments in this direction may again rest with the complicated 
form the equations of motion take if one starts from the idea of defining the 
potential in terms of the intrinsic distance in the framework of differential 
geometry. Such complications might have discouraged all the attempts to 
generalize the problem to more than two bodies. 

1.3. Our approach. The present paper overcomes the above-mentioned diffi- 
culties encountered in defining a meaningful n-body problem prone to the same 
mathematical depth achieved in the classical case, by replacing the differential- 
geometric approach used for n = 2 in the case of the cotangent potential with 
the variational method of constrained Lagrangian dynamics. Also, the techni- 
cal complications that arise in understanding the motion within the standard 



A vivid description of the collaboration between Einstein and Infeld appears in |37| . 
^One of us (Erensto Perez-Chavela) , together with his student Luis Franco-Perez, recently 
analyzed a restricted 3-body problem in S^, |31| . in a more restrained context than the one 
we provide here. 
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models of the Bolyai-Lobachevsky plane (the Klein-Beltrami disk, the Poincare 
upper-half-plane, and the Poincare disk) are bypassed through the less known 
Weierstrass hyperboloidal model (see Appendix), which often provides analo- 
gies with the results we obtain in the spherical case. This model also allows 
us to use hyperbolic rotations — a class of isometries — to put into the evidence 
some unexpected solutions of the equations of motion. 

The history of the problem shows that there is no unique way of extending 
the classical idea of gravitation to spaces of constant curvature, but that the 
cotangent potential is the most natural candidate. Therefore we take this 
potential as a starting point of our approach, though some of our results — as 
for example Saari's conjecture in the geodesic case — do not use this potential 
explicitly, only its property of being a homogenous function of degree zero. 

Our generalization recovers the Newtonian law when the curvature is zero. 
Moreover, it provides a unified context, in which the potential varies continu- 
ously with the curvature k. The same continuity occurs for the basic results 
when the curvature tends to zero. For instance, the set of closed orbits of the 
Kepler problem on non-zero-curvature surfaces tends to the set of ellipses in 
the Euchdean plane when /t — > (see, e.g., [H] or [Sj). 

2. Summary of results 

2.1. Equations of motion. In Section 3, we extend the Newtonian poten- 
tial of the n-body problem to spaces of constant curvature, k, for any finite 
dimension. For k 7^ 0, the potential turns out to be a homogeneous function 
of degree zero. We also show the existence of an energy integral as well as of 
the integrals of the angular momentum. Like in general relativity, there are no 
integrals of the center of mass and linear momentum. But unlike in relativity, 
where — in the passage from continuous matter to discrete bodies — the fact 
that forces don't cancel at the center of mass leads to difficulties in defining 
infinitesimal sizes for finite masses, [Sj, we do not encounter such problems 
here. We assume that the laws of classical mechanics hold for point masses 
moving on manifolds, so we can apply the results of constrained Lagrangian 
dynamics to derive the equations of motion. Thus two kinds of forces act 
on bodies: (i) those given by the mutual interaction between particles, repre- 
sented by the gradient of the potential, and (ii) those that occur due to the 
constraints, which involve both position and velocity terms. 

2.2. Singularities. In Section 4 we focus on singularities, and distinguish 
between singularities of the equations of motion and solution singularities. 
For any k 7^ 0, the equations of motion become singular at collisions, the 
same as in the Euclidean case. The case k > 0, however, introduces some new 
singularities, which we call antipodal because they occur when two bodies are 
at the opposite ends of a diameter of the sphere. 
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The set of singularities is endowed with a natural dynamical structure. 
When the motion of three bodies takes place along a geodesic, solutions close 
to binary collisions and away from antipodal singularities end up in collision, 
so binary collisions are attractive. But antipodal singularities are repulsive in 
the sense that no matter how close two bodies are to an antipodal singularity, 
they never reach it if the third body is far from a collision with any of them. 

Solution singularities arise naturally from the question of existence and 
uniqueness of initial value problems. For nonsingular initial conditions, stan- 
dard results of the theory of differential equations ensure local existence and 
uniqueness of an analytic solution defined in some interval [0,t"'"). This solu- 
tion can be analytically extended to an interval [0,t*), with < t^ < t* < oo. 
If t* = CX3, the solution is globally defined. If t* < oo, the solution is called 
singular and is said to have a singularity at time t*. 

While the existence of solutions ending in collisions is obvious for any value 
of K, the occurrence of other singularities is not easy to demonstrate. Nev- 
ertheless, we prove that some hybrid singular solutions exist in the 3-body 
problem with n > 0. These orbits end up in finite time in a collision-antipodal 
singularity. Whether other types of non-collision singularities exist, like the 
pseudocollisions of the Euclidean case, remains an open question. The main 
reason why this problem is not easy to answer rests with the nonexistence of 
the center-of-mass integrals. 

Another class of solutions connected to collision-antipodal configurations 
is particularly interesting. We show that, for n = 3, there are orbits that 
reach such a configuration at some instant t* but remain analytic at this point 
because the forces and the velocities involved remain finite at t*. Such a motion 
can, of course, be analytically continued beyond t*. This is the first example 
of a natural non-singularity collision. 

2.3. Relative equilibria. The rest of this paper, except for the Appendix, 
focuses on the results we obtained in S^ and H^, mainly because these two 
surfaces are representative for the cases k > and k < 0, respectively. Indeed, 
the results we proved for these surfaces can be extended to different curvatures 
of the same sign by a mere change of factor. 

Sections 5 and 6 deal with relative equilibria in S^ and H^. In S^ we only 
have elliptic relative equilibria. Instead, the relative equilibria in H^ are of 
two kinds: elliptic relative equilibria, generated by elliptic rotations, and hy- 
perbolic relative equilibria, generated by hyperbolic rotations (see Appendix). 
Parabolic relative equilibria, generated by parabolic rotations, do not exist. 

Some of the results we obtain in S^ have analogues in H^; others are specific 
to each case. Theorems [6] and [TOl for instance, are dual to each other, whereas 
Theorem |2] takes place only in S^. The latter identifies a class of fixed points 
of the equations of motion. More precisely, we prove that if an odd number 
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n of equal masses are placed, initially at rest, at the vertices of a regular 
ra-gon inscribed in a great circle, then the bodies won't move. The same 
is true for four equal masses placed at the vertices of a regular tetrahedron 
inscribed in S^, but — due to the occurrence of antipodal singularities — fails to 
hold for the other regular polyhedra: octahedron (6 bodies), cube (8 bodies), 
dodecahedron (12 bodies), and icosahedron (20 bodies), as well as in the case 
of geodesic n-gons with an even number of bodies. 

Theorem O shows that there are no fixed points for n bodies within any 
hemisphere of S^. Its hyperbolic analogue, stated in Theorem [HI proves the 
nonexistence of fixed points in H^. These two results are in agreement with 
the Euclidean case in the sense that the n-body problem has no fixed points 
within distances, say, not larger than the ray of the visible universe. 

It is also natural to ask whether fixed points can generate relative equilibria. 
Theorem H] shows that if n masses mi,m2, . . . , m„ lie initially on a great circle 
of S^ such that the mutual forces are in equilibrium, then any uniform rotation 
applied to the system generates a relative equilibrium. 

Theorem [S] states that the only way to generate an elliptic relative equilib- 
rium from an initial n-gon configuration taken on a great circle, as in Theorem 
[21 is to assign suitable velocities in the plane of the n-gon. So a regular polygon 
of this kind can rotate only in a plane orthogonal to the rotation axis. 

Theorem [6] and its hyperbolic analogue, Theorem [TOl show that n-gons of 
any admissible size can rotate on the same circle, both in S^ and H^. Again, 
these results agree with the Euclidean case. But something interesting happens 
with the equilateral (Lagrangian) solutions. Unlike in Euclidean space, elliptic 
relative equilibria moving in the same plane of R^ can be generated only when 
the masses move on the same circle and are therefore equal, as we prove 
in Theorems [7] and [TTJ Thus Lagrangian solutions with unequal masses are 
specific to the Euclidean case. 

Theorems [8] and [12] show that analogues to the coUinear (Eulerian) orbits 
in the 3-body problem of the classical case exist in S^ and H^, respectively. 
While nothing surprising happens in H^, where we prove the existence of such 
solutions of any size, an interesting phenomenon takes place in S^. Assume 
that one body lies on the rotation axis (which contains one height of the tri- 
angle), while the other two are at the opposite ends of a rotating diameter on 
some non-geodesic circle of S^. Then elliptic relative equilibria exist while the 
bodies are at initial positions within the same hemisphere. When the rotat- 
ing bodies are placed on the equator, however, they encounter an antipodal 
singularity. Below the equator, solutions exist again until the bodies form an 
equilateral triangle. By Theorem [Sj any n-gon with an odd number of sides 
can rotate only in its own plane, so the (vertical) equilateral triangle is a fixed 
point but cannot lead to an elliptic relative equilibrium. If the rotating bodies 
are then placed below the equilateral position, solutions fail to exist. But the 
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masses don't have to be all equal. Eulerian solutions exist if, say, the non- 
rotating body has mass m and the other two have mass M. If M > Am, these 
orbits occur for all ;i; 7^ 0. Again, these results prove that, as long as we do 
not exceed reasonable distances, such as the ray of the visible universe, the 
behavior of elliptic relative equilibria lying on a rotating geodesic is similar to 
the one of Eulerian solutions in the Euclidean case. 

We then study hyperbolic relative equilibria around a point and along a 
(usually non-geodesic) hyperbola. Theorem [12] proves that such orbits do not 
exist on fixed geodesies of H^, so the bodies cannot chase each other along a 
geodesic while maintaining the same initial distances. But Theorem UM proves 
the existence of hyperbolic relative equilibria in H^ for three equal masses. The 
bodies move along hyperbolas of the hyperboloid that models H^ remaining 
all the time on a moving geodesic and maintaining the initial distances among 
themselves. These orbits rather resemble fighter planes flying in formation 
than celestial bodies moving under the action of gravity alone. The result also 
holds if the mass in the middle differs from the other two. The last result of 
this section. Theorem \T5[ shows that parabolic relative equilibria do not exist. 

2.4. Saari's conjecture. Our extension of the Newtonian n-body problem 
to spaces of constant curvature also reveals new aspects of Saari's conjecture. 
Proposed in 1970 by Don Saari in the Euclidean case, Saari's conjecture claims 
that solutions with constant moment of inertia are relative equilibria. This 
problem generated a lot of interest from the very beginning, but also several 
failed attempts to prove it. The discovery of the figure eight solution, which 
has an almost constant moment of inertia, and whose existence was proved in 
2000 by Alain Chenciner and Richard Montgomery, [2], renewed the interest 
in this conjecture. Several results showed up not long thereafter. The case 
n = 3 was solved in 2005 by Rick Moeckel, [56]; the coUinear case, for any 
number of bodies and the more general potentials that involve only mutual 
distances, was settled the same year by the authors of this paper, [25]. Saari's 
conjecture is also connected to the Chazy-Wintner-Smale conjecture, [70], [76] , 
which asks to determine whether the number of central configurations is finite 
for n given bodies in Euclidean space. 

Since relative equilibria have elliptic and hyperbolic versions in H^, Saari's 
conjecture raises new questions for k < 0. We answered them in Theorem [TBI of 
Section 7, when the bodies are restrained to a geodesic that rotates elliptically 
or hyperbolically. 

An Appendix in which we present some basic facts about the Weierstrass 
model of the hyperbolic plane, together with some historical remarks, closes 
our paper. We suggest that readers unfamiliar with this model take a look at 
the Appendix before getting into the technical details related to our results. 
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2.5. Some physical remarks. Does our gravitational model have any con- 
nection with the physical reality? Since there is no unique extension of the 
Newtonian n-body problem to spaces of constant curvature, is our generaliza- 
tion meaningful from the physical point of view or does it lead only to some 
interesting mathematical properties? 

We followed the tradition of the cotangent potential, which seems the most 
natural candidate. But since the debate on the nature of the physical space is 
open, the only way to justify this model is through mathematical results. As 
we will further argue, not only that the properties we obtained match the Eu- 
clidean ones, but they also provide a classical explanation of the cosmological 
scenario, in agreement with the basic conclusions of general relativity. 

But before getting into the physical aspect, let us remark that our model 
is based on mathematical principles, which lead to a meaningful physical in- 
terpretation. As we already mentioned, the cotangent potential preserves two 
fundamental properties: (i) it is harmonic for the one-body problem and (ii) 
it generates a central field in which all bounded orbits are closed. Other re- 
sults that support the cotangent potential are based on the idea of central 
(or gnomonic) projection, p. By taking the central projection on the sphere 
for the planar Kepler problem, Paul Appell obtained the cotangent potential. 
This idea can be generalized by projecting the planar Kepler problem to any 
surface of revolution, as one of us (Manuele Santoprte) proved, [6T] . 

In 1992, Kozlov and Harin showed that the only central potential that sat- 
isfies the fundamental properties (i) and (ii) in S^ and has meaning in celestial 
mechanics is the cotangent of the distance, [13]. This fact had been known 
to Infeld for the quantum mechanical version of the potential, [36] • But since 
any continuously differentiable and non-constant harmonic function attains no 
maximum or minimum on the sphere, the existence of two distinct singularities 
(the collisional and the antipodal — in our case) is not unexpected. And though 
a force that becomes infinite for points at opposite poles may seem counterin- 
tuitive in a gravitational framework, it explains the cosmological scenario. 

Indeed, while there is no doubt that n point masses ejecting from a total 
collapse would move forever in spaces with k < for large initial conditions, 
in agreement with general relativity, it is not clear what happens for k > 
0. But the energy relation 0221) shows that, in spherical space, the current 
expansion of the universe cannot last forever. For a fixed energy constant, 
h, the potential energy, — [/, would become positive and very large if one or 
more pairs of particles were to come close to antipodal singularities. Therefore 
in a homogeneous universe, highly populated with non-colliding particles, the 
system could never expand beyond the equator (assuming that the initial 
ejection took place at one pole). No matter how large (but fixed) the energy 
constant is, when the potential energy reaches the value /i, the kinetic energy 
becomes zero, so all the particles stop simultaneously and the motion reverses. 
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Thus, for K > 0, the cotangent potential recovers the growth of the system to 
a maximum size and the reversal of the expansion independently on the value 
of the energy constant. Without antipodal singularities, the reversal could take 
place only for certain initial conditions. This conclusion is reached without 
introducing a cosmological force and differently from how it was obtained in 
the classical model proposed by Elie Cartan, |Tl], [12], and shown by Frank 
Tipler to be as rigorous as Friedmann's cosmology, [72], [73] . 

Another result that suggests the validity of the cotangent potential is the 
nonexistence of fixed points. They don't show up in the Euclidean case, and 
neither do they appear in this model within the observable universe. The 
properties we proved for relative equilibria are also in agreement with the 
classical n-body problem, the only exception being the Lagrangian solutions 
for K ^ 0, which, unlike in the Euclidean case, must have equal masses and 
move on the same circle. This distinction adds to the strength of the model 
because, even in the Euclidean case, the arbitrariness of the Lagrangian so- 
lutions is a peculiar property. At least two arguments support this point of 
view. First, relative equilibria generated from regular polygons, except the 
equilateral triangle, exist only if the masses are equal. The second argument 
is related to central configurations, which generate relative equilibria in the 
Euclidean case. One of us (Florin Diacu) proved that among attraction forces 
given by symmetric laws of masses, 7(mj,mj) = '~f{mj,mi), equilateral cen- 
tral configurations with unequal masses occur only when '~f{mi,'mj) = cmimj, 
where c is a nonzero constant, [20]. Since for n ^ relative equilibria are 
equilateral only if the masses are equal means that Lagrangian solutions of 
arbitrary masses characterize the Euclidean space. 

Such orbits exist in nature, the best known example being the equilateral 
triangle formed by the Sun, Jupiter, and the Trojan asteroids. Therefore our 
result reinforces the fact that space is Euclidean within distances comparable 
to those of our solar system. This fact was not known during the time of 
Gauss, who apparently tried to determine the nature of space by measuring 
the angles of triangles having the vertices some tens of kilometers apart |j 
Since we cannot measure the angles of cosmic triangles, our result opens up 
a new possibility. Any evidence of a Lagrangian solution involving galaxies 
(or clusters of galaxies) of unequal masses, could be used as an argument for 
the flatness of the physical space for distances comparable to the size of the 
triangle. Similarly, hyperbolic relative equilibria would show that space has 
negative curvature. 



Arthur Miller argues that these experiments never took place, 
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3. Equations of motion 

We derive in this section a Newtonian n-body problem on surfaces of con- 
stant curvature. The equations of motion we obtain are simple enough to allow 
an analytic approach. At the end, we provide a straightforward generalization 
of these equations to spaces of constant curvature of any finite dimension. 

3.1. Unified trigonometry. Let us first consider what, following [5], we will 
call trigonometric K-functions, which unify elliptical and hyperbolic trigonom- 
etry. We define the K-sine, sn^, as 

{fi;~^/^sinK^/^x if k > 
X if K = 
{-k)-^/^ smh{- Ky/^x if ft<0, 

the K-cosine, csn^, as 

{cosK^/^x if K > 
1 if K = 
cosh(-K)^/^x if K < 0, 

as well as the K-tangent, tn^, and k- cotangent, ctn^, as 

tn^[x) := -— and ctn^^x) := -— , 

csn«;(x) sn«;(a;) 

respectively. The entire trigonometry can be rewritten in this unified context, 
but the only identity we will further need is the fundamental formula 

K sn^(x) + csn^(x) = 1. 

3.2. Differential-geometric approach. In a 2-dimensional Riemann space, 
we can define geodesic polar coordinates, (r, </>), by fixing an origin and an 
oriented geodesic through it. If the space has constant curvature n, the range 
of r depends on k; namely r G [0, 7r/(2fi:^'^)] for k, > and r G [0, oo) for 
K < 0; in all cases, (j) G [0, 27i]. The line element is given by 

dsl = dr^ + snl{r)d(j)^. 

In S^,R^, and H^, the line element corresponds to k = 1,0, and —1, respec- 
tively, and reduces therefore to 

dsi = dr"^ + (sin^ r)(i<^^, dsQ = dr"^ + r'^dcf?^ and ds'^_i = dr"^ + (sinh r)d(l)'^. 

In [9] , the Lagrangian of the Kepler problem is defined as 

where Vr and Vfp represent the polar components of the velocity, and —U is the 
potential, where 

Ujr) = G ctnjr) 
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is the force function, G > being the gravitational constant. This means that 
the corresponding force functions in S^,R^, and H^ are, respectively, 

Ui{r) = Gcotr, Uo{r) = Gr^^, and f/-_i(r) = Gcothr. 

In this setting, the case k = separates the potentials with k > and k < 
into classes exhibiting different qualitative behavior. The passage from k > 
to K, < through k, = takes place continuously. Moreover, the potential 
is spherically symmetric and satisfies Gauss's law in a 3-dimensional space 
of constant curvature k. This law asks that the flux of the radial force field 
across a sphere of radius r is a constant independent of r. Since the area of 
the sphere is 47rsn^(r), the flux is 47rsn^(r) x ■^[/^{r), so the potential satisfies 
Gauss's law. As in the Euclidean case, this generalized potential does not 
satisfy Gauss's law in the 2-dimensional space. The results obtained in [9] 
show that the force function U^ leads to the expected conic orbits on surfaces 
of constant curvature, and thus justify this extension of the Kepler problem 
to K ^ 0. 

3.3. The potential. To generalize the above setting of the Kepler problem 
to the 77,-body problem on surfaces of constant curvature, let us start with 
some notations. Consider n bodies of masses mi, . . . , m„ moving on a surface 
of constant curvature k. When k > 0, the surfaces are spheres of radii k~^/^ 
given by the equation x'^ + y'^ + z'^ = k,~^; ioi k, = 0, we recover the Euclidean 
plane; and if k < 0, we consider the Weierstrass model of hyperbolic geometry 
(see Appendix) , which is devised on the sheets with z > of the hyperboloids 
of two sheets x'^ + y'^ — z"^ = k,~^. The coordinates of the body of mass rrii are 
given by Qj = {xi,yi,Zi) and a constraint, depending on k, that restricts the 
motion of this body to one of the above described surfaces. 
In this paper, Vq- denotes either of the gradient operators 

Vq, = {d^^,dy^,d,J, for k > 0, or Vq, = {d^^,dy.,-d,J, for k < 0, 

with respect to the vector q^, and V stands for the operator (Vq^, . . . , Vq,J. 
For a = {ttx, tty, ttz) and b = {bx, by, b^) in R"^, we define a b as either of the 
inner products 

a ■ b := {a^bx + CLyby + ttzb^) for k > 0, 

a H b := {ttxbx + ayby — ttzbz) for k < 0, 

the latter being the Lorentz inner product (see Appendix). We also define 
a ® b as either of the cross products 

a X b := {aybz - azby, azbx - axbz, axby - aybx) for k > 0, 
a K b := {ttybz - azby, azbx - axbz, aybx - axby) for k <0. 
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The distance between a and b on the surface of constant curvature n is then 
given by 

fi;~^/^cos~^(fi:a- b), k>0 

d^isL, b) := ■^ |a — b|, k = 

where the vertical bars denote the standard Euchdean norm. In particular, 
the distances in S^ and H^ are 

di{a, b) = cos""'^(a ■ b), (i_i(a, b) = cosh~^(— a El b), 

respectively. Notice that do is the limiting case of d^ when k ^ 0. Indeed, 
for both K, > and k < 0, the vectors a and b tend to infinity and become 
parallel, while the surfaces tend to an Euclidean plane, therefore the length of 
the arc between the vectors tends to the Euclidean distance. 

We will further define a potential in R^ if k > 0, and in the 3-dimensional 
Minkowski space A1^ (see Appendix) if k < 0, such that we can use a varia- 
tional method to derive the equations of motion. For this purpose we need to 
extend the distance to these spaces. We do this by redefining the distance as 

■^-i/^cos-i^^s:^ ;.>0 

d^,{a,h) := { \a-h\, k = 

-)"''' -^h-7^sS?=bH5''^<0- 

Notice that this new definition is identical with the previous one when we 
restrict the vectors a and b to the spheres x^+y'^+z^ = k~^ or the hyperboloids 
x"^ + y"^ — z'^ = K~^, but is also valid for any vectors a and b in R'^ and Ai^, 
respectively. 

From now on we will rescale the units such that the gravitational constant 
G is 1. We thus define the potential of the n-body problem as the function 
-f/^(q), where 

Uf^ifi) ■= -"^ ^ mimjCtn«,(4(qi,qj)) 

i=l j=l,JT^i 

stands for the force function, and q = (qi,...,q„) is the configuration of 
the system. Notice that ctno((io(qi, qj)) = |qi — qj|~^, which means that we 
recover the Newtonian potential in the Euclidean case. Therefore the potential 
Uk varies continuously with the curvature k. 

Now that we defined a potential that satisfies the basic continuity condition 
we required of any extension of the n-body problem beyond the Euclidean 
space, we will focus on the case k 7^ 0. A straightforward computation shows 
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that 

-. n n UliTfljiCrKj ' - — , 

(1) [/4q) = -2^ ^ ,2 "' ^^0' 



2 

1=1 j=l,jj^i 



a-a\ ^^^^^^ 



VKqiOqi-^KqjOqj 



where 



-1, for K > 
-1, for K < 0. 



3.4. Euler's formula. Notice that U^ifjci) = U^i^q) = rpUi^{(\} for any rj ^ 0, 
which means that the potential is a homogeneous function of degree zero. But 
for q in R^", homogeneous functions F : R^" — > R of degree a satisfy Euler's 
formula, q ■ V-F(q) = aF{q). With our notations, Euler's formula can be 
written as q VF(q) = aF(q). Since a = for 11^ with k 7^ 0, we conclude 
that 

(2) q0W«(q) = O. 

We can also write the force function as Uf^{q) = \ Ym=i ^kIq*); where 

TTii \ sr^ V^qiOqiV'ijQ'ij 



, v^req^Oqly^iq/Dqj, 

are also homogeneous functions of degree 0. Applying Euler's formula for 
functions F : R^ — > R, we obtain that q^ © Vq.?7*(q) = 0. Then using the 
identity Vq-f/K(q) = Vq-f/*(qj), we can conclude that 

(3) qi0Vq,f/«(q) = O, z = l,...,n. 

3.5. Derivation of the equations of motion. To obtain the equations of 
motion for k 7^ 0, we will use a variational method applied to the force function 
([1]). The Lagrangian of the n-body system has the form 

L,(q,q)=T«(q,q) + f/,(q), 

where TK(q, q) := \Yl'i=i^i{^i © qj)('^qi © Qj) is the kinetic energy of the 
system. (The reason for introducing the factors /tqjOqj = 1 into the definition 
of the kinetic energy will become clear in Section 3.8.) Then, according to the 
theory of constrained Lagrangian dynamics (see, e.g., [32]), the equations of 
motion are 
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where /^ = qjOq^ — k^^ is the function that gives the constraint /^ = 0, which 
keeps the body of mass rrii on the surface of constant curvature k, and A];, is 
the Lagrange multipher corresponding to the same body. Since q^ q^ = k"^ 
imphes that qj qj = 0, it follows that 

d I dLf^ \ . 

-t; I -^ I = miqi[KCii qj) + 2mi[K<ii qi) = mjqi. 

This relation, together with 

dL ~ 

Cqi 
implies that equations (jll) are equivalent to 

(5) uiiiii - rriiK^qi qi)qi - Vci,U^,{q) - 2X'^{t)qi = 0, i = l,...,n. 
To determine AJ,, notice that = /^ = 2qj q^ + 2(qj qj), so 

(6) qi0qi = -qi0qi- 

Let us also remark that 0-multiplying equations ([5]) by q^ and using ([3]), we 
obtain that 

mi{qi qi) - mjfi;(qi q^) - q^ Vq,f/K(q) = 2A*^qi q^ = 2fi:"^A*^, 

which, via IQ, implies that A^ = — Kmj(qj q,). Substituting these values of 
the Lagrange multipliers into equations ([5]), the equations of motion and their 
constraints become 

(7) m^qi = Vq^f/^(q) -miK(qi0qi)qi, q^ q^ = fi:"\ k 7^ 0, 

i = 1, . . . ,n. 
The qj-gradient of the force function, obtained from ([1]), has the form 

, si /if K^qiOqj 



Vq,f/,(q) = 2^ ^ — ^, /€ ^ 0, 



i=ijVi 



o-_o- . ^5u0q2_ 



VftqiOqi-y/KqjOqj 

and using the fact that /tqj q,j = 1, we can write this gradient as 



mjmj^anf/'^ [q,- - (/tq^ q^ - ) q^' 



(9) v,f/.(q) = E "^-"^^-^^"^ iq.-i^q^Qq,-jq.j ^ ^ _^ ^^ 



Sometimes we can use the simpler form ([9]) of the gradient, but whenever we 
need to exploit the homogeneity of the gradient or have to differentiate it, 
we must use its original form ([8]). Thus equations ([7]) and ([8]) describe the 
n-body problem on surfaces of constant curvature for k 7^ 0. Though more 
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complicated than the equations of motion Newton derived for the Euchdean 
space, system ([7]) is simple enough to allow an analytic approach. Let us first 
provide some of its basic properties. 

3.6. First integrals. The equations of motion have the energy integral 

(10) T«(q,p)-t/«(q) = /i, 

where, recall, Tf,(q, p) := | Yl^=i '^j~^(Pi © Pi){f^^i © ^.i) is the kinetic energy, 
p := (pi, . . . , p„) denotes the momentum of the n-body system, with pj : = 
miC[i representing the momentum of the body of mass rrii, i = 1, . . . ,n, and h 
is a real constant. Indeed, ©-multiplying equations ([7]) by qj, we obtain 

^m^qi © qi = [Vq,U,,{q)] © q^ - ^miK(qi © qi)qi © q^ = —U,,{q{t)). 

i=l i=l 

Then equation ( ITOi) follows by integrating the first and last term in the above 
equation. 

The equations of motion also have the integrals of the angular momentum, 

n 

(11) ^qi©pi = c, 

where c is a constant vector. Relations (TTTj) follow by integrating the identity 
formed by the first and last term of the equations 

^ V^ V^ mimj-((TK)3/2q. ®q. 

(12) > niiCii © qi = ) ) '— -—j- 



-E 



Emimj[aKfl^[Ki\i © q^) / • ^ • 
'— —r^ - m^K q, © q^ 



qj © q^ = 0, 



obtained if ©-multiplying the equations of motion ([7]) by qj. The last of 
the above identities follows from the skew-symmetry of © and the fact that 

qj ©qi = 0, z = 1,. . . ,n. 

3.7. Motion of a free body. A consequence of the integrals of motion is 
the analogue of the well known result from the Euclidean space related to 
the motion of a single body in the absence of any gravitational interactions. 
Though simple, the proof of this property is not as trivial as in the classical 
case. 

Proposition 1. A free body on a surface of constant curvature is either at 
rest or it moves uniformly along a geodesic. Moreover, for k > 0, every orbit 
is closed. 
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Proof. Since there are no gravitational interactions, the equations of motion 
take the form 

(13) q = -fi;(q©q)q, 

where q = (x, y, z) is the vector describing the position of the body of mass 
m. If q(0) = 0, then q(0) = 0, so no force acts on m. Therefore the body will 
be at rest. 

If q(0) 7^ 0, q(0) and q(0) are coUinear, having the same sense if k < 0, but 
the opposite sense if k > 0. So the sum between q(0) and q(0) pulls the body 
along the geodesic corresponding to the direction of these vectors. 

We still need to show that the motion is uniform. This fact follows obviously 
from the integral of energy, but we can also derive it from the integrals of the 
angular momentum. Indeed, for k > 0, these integrals lead us to 

c = (q X q) ■ (q X q) = (q ■ q)(q ■ q) sin^ a, 

where c is the length of the angular momentum vector and a is the angle 
between q and q (namely vr/2). So since q ■ q = k~^, we can draw the 
conclusion that the speed of the body is constant. 
For K < 0, we can write that 



c = (q Kl q) H (q K q) 



K-i 
qHq 



-K ^qElq. 



qHq qHq 
qHq qHq 

Therefore the speed is constant in this case too, so the motion is uniform. Since 
for K > the body moves on geodesies of a sphere, every orbit is closed. H 




3.8. Hamiltonian form. The equations of motion ([7]) are Hamiltonian. In- 
deed, the Hamiltonian function H,^ is given by 

= \ Er=i ^'i^ipi Pi)('«qi qi) - u^{(\), 

K~^, K 7^ 0, i = 1, . . . ,n. 

Equations ([5]) thus take the form of a 6n-dimensional first order system of 
differential equations with 2n constraints, 

{q* = Vp^H^{q,p) = m'^pi, 
Pi = -Vq,i/.(q, p) = Vq,f/«(q) - m-^K{pi pi)qi, 
qi © qi = K-\ qj © Pi = 0, k 7^ 0, z = 1, . . . , n. 

It is interesting to note that, independently of whether the kinetic energy is 
defined as 

-, n 1 " 

TM ■= 2^^i^'P' ® Pi °^ ^^(q.p) := 2 5Z"^r^(Pi0pO(«^qi0qi)) 

i=l i=l 

(which, though identical since Kq^ © q, = 1, does not come to the same thing 
when differentiating T^), the form of equations ([7]) remains the same. But 
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in the former case, system (ITj) cannot be put in Hamiltonian form in spite of 
having an energy integral, while in the former case it can. This is why we 
chose the latter definition of T^. 

These equations describe the motion of the n-body system for any k 7^ 0, 
the case k = corresponding to the classical Newtonian equations. The 
representative non- zero-curvature cases, however, are k = 1 and k = —1, 
which characterize the motion for k > and k < 0, respectively. Therefore 
we will further focus on the n-body problem in S^ and H^. 



3.9. Equations of motion in S^. In this case, the force function ([T]) takes 
the form 



(15) 






*=i j='i^,j¥=i , I 



qi-qj 



.vqT^v/qrqJ, 

while the equations of motion ([7]) and their constraints become 
(16) miiii = Vq,t/i(q)-mi(qj-qj)qi, qj-qj = l, qi-qi = 0, i = 1, 



,n. 



In terms of coordinates, the equations of motion and their constraints can 
be written as 



(17) 



rriiyi 
rriiZi 



^-m,{x^ + yf + zf)x„ 
^ ~ m,{xj + yf + zf)yu 
^-m,{x^ + yf + zf)zi, 



{xj + yf + zf 



1, XiX'i ~r yiyi \ Z^Zi 

and by computing the gradients they become 



0, i 



l,...,n. 



Xi 



En 



J^i 



x^Xj+y^yj+ZiZj 



V^f+yf+^h/^J+yf+- 









3 "] 3 



3/2 



[xi + yf + zt)xi 



En 



a/"?+^' + ^V"j+'^j + 



1- 






Tj2 



- [Xt 



yf 



zf)yu 



X- 



En 



yf + zf = 1, 



yM+yf+^h/^J+yl 



XiXj+yiVj+ZiZj 



\A!+^!+^ 



y/x^+y^ + z^ 



-1 3/2 



J J J 

XiXi + yiyi ~r ZiZi U 



[il 



1,- 



2/' + 



,n. 



zr Zi 
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Since we will neither need the homogeneity of the gradient, nor we will we 
differentiate it, we can use the constraints to write the above system as 



(19) 



™ _ V" mj[xj-{xiXj+yiyj+ZiZj)xi] _ / ■ 2 , -2 , •2\™ 

.. _ Y^n rrij [yj-{xiXj+yiyj+ZiZj)yi] _ / . 2 -2 , ■2\.. 



mj[zj-{xiXj+yiyj+ZiZj)zi] _ / . 2 , -2 , -2 
jVi [l-{x,Xj+yiyj+z,z,)2]y2 ^-^i ^ ^i ^ ^i 

-^-1, XiXi + yiyi + ZiZi = Q, z = l,.. 
The Hamiltonian form of the equations of motion is 



Ixf + yf + zt 



n. 



(20) 



"^i Pi, 

m.imj [qj-(qt-q.,)qt] _ 
^3 = ^,3+i [l-(qi-q,-)2]3/' 

qi ■ Qj = 1, Qi ■ Pi = 0, K 7^ 0, 



En 
7 = 1, 



-1/ 



2 = 1, . . . ,n. 
Consequently the integral of energy has the form 



(21) 






rriim 



<\i<ij 



'-' vqi-qvij'ij 



i=l 



2/i, 



1- 



q»qj 



V^rqiv^q^qJ 



which, via q^ ■ q^ 

n 

(22) E 



1, -i = 1, . . . , n, becomes 



?Ti, 



"1, 



Pi ■ Pi 



EE 



mirrijfii ■ q^- 



and the integrals of the angular momentum take the form 



2/i, 



(23) 



E qi X Pi = c. 



i=l 



Notice that sometimes we can use the simpler form (I22p of the energy integral, 
but whenever we need to exploit the homogeneity of the potential or have to 
differentiate it, we must use the more complicated form (1211) . 



3.10. Equations of motion in H^. In this case, the force function ([T]) takes 
the form 

„ n mm- — ^ 

(24) t^-i(q) = -^EE 



■rriimj 



^-q, Hq, -y/-qj Sqj 



i=i j=ijyi 



qiQqj 



^^-qiSqi^-qjQqj 

SO the equations of motion and their constraints become 

(25) miCJi = Vq,,f/_i(q) + mi(qi B qi)qi, qi H qi = -1, qi B qi 



n. 



The n-Body Problem in Spaces of Constant Curvature 



21 



In terms of coordinates, the equations of motion and their constraints can 
be written as 



(26) 



rriiiji = 
rriiZi = ■ 

{x^ + yf 






= -1 



zf)yi 



yf - zf)zi 



z? 



•^i-^i ~r ytlji ZiZi U, % 



1,- 



and by computing the gradients they become 



Xj^Xj +y{yj ^i^j 



En 



v^- 



2 /_^2_ 2 + ^2 

i V J J J 



x^Xj-^-y^yj ■^i^j 



. W""=»"''?+'^?A/"''i"^i+' 



^i^j+ViVj 



3/2 



+ (i? + iji - zt)xi, 






(27) { 



j-j. - V" - 



V""?~^? + ''?A/-"h^i+". 



Xj^Xj -f-yiyj ^i^j 



z- = V 



^-xj-y^+.2^_^2_y2+^ 



TJ72 + (xl + y»^ - zf)yi 



m,- 



'j='i-,j¥=i 



A + 2/- - 2;2 



^,,2_^2 + ,2^_,2_^2 + , 



Xj^Xj -j-yj^yj Zj^Zj 






3/2 



+ {it + fi - zt)zi, 



,n. 



For the same reasons described in the previous subsection, we can use the 
constraints to write from now on the above system as 



(28) 



Xi 



mj[xj + {xiXj+yiyj-ZiZj)xi\ / -2 , -2 _ •2\™ 

V" mj [yj+ixiXj+yiyj-ZiZj)yi] / -2 , -2 _ •2\„, 

Z^7=l,i9^j [('7....a-..-l-„.-7,..-^.-^..12_ii3/2 "T V-^j -T yj ^i ;(/i, 



En 
.7 = 1,. 



d='i-,j¥=i [(xiXj+yiyj-ZiZj)2-l]3/2 
mj[zj+{xiXj+yiyj-ZiZj)zi] 



En 



^i [{xiXj+y^yj- 



ZiZiY-\\^l^ 



A^yi 



■1, XiXi -]- yiyi 



^i ^i 



{il + yl 

= 0, i = 



l,...,n. 



The Hamiltonian form of the equations of motion is 



% 



"^i Pi 



(29) 



_ sr^n m^ m j [qj + {eg Eqj ) q^ 

Pi - Z^j = ij^j [{q,mqj)2-l]3/2 



qjHqi 



-1, qiE]pi = 0, 



m. 



(pi H pi)qj 



l,...,n. 
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Consequently the integral of energy takes the form 

n n n TfliTflj — , = — , = 

(30) ^m-'(p.0p.) + 5: Z -^^S^^ = 2h. 

which, via qj H qj = — 1, z = 1, . . . , ra, becomes 



1 



n n 



(31) E'".-'(P.0P.) + E E ^pm^^2K 
and the integrals of the angular momentum can be written as 

n 

(32) ^q,Kpi = c. 

Notice that sometimes we can use the simpler form (131 p of the energy integral, 
but whenever we need to exploit the homogeneity of the potential or have to 
differentiate it, we must use the more complicated form ( l30l) . 

3.11. Equations of motion in S'^ and H^. The formalism we adopted in 
this paper allows a straightforward generalization of the n-body problem to 
S^ and H'^ for any integer /i > 1. The equations of motion in /^-dimensional 
spaces of constant curvature have the form ([7]) for vectors qj and q^- of R^+^ 
constrained to the corresponding manifold. It is then easy to see from any 
coordinate-form of the system that S"^ and H'^ are invariant sets for the equa- 
tions of motion in S'* and H^, respectively, for any integer u < fi. 

Indeed, this is the case, say, for equations flT9|) . if we take Xi{0) = 0, Xi(0) = 
0, i = 1, ... ,77,. Then the equations of Xi are identically satisfied, and the 
motion takes place on the circle y"^ + z"^ = 1. The generalization of this idea 
from one component to any number u of components in a (/i + l)-dimensional 
space, with u < fi, is straightforward. Therefore the study of the n-body 
problem on surfaces of constant curvature is fully justified. 

The only aspect of this generalization that is not obvious from our formalism 
is how to extend the cross product to higher dimensions. But this extension 
can be done as in general relativity with the help of the exterior product. 
However, we will not get into higher dimensions in this paper. Our further 
goal is to study the 2-dimensional case. 

4. Singularities 

Singularities have always been a rich source of research in the theory of 
differential equations. The n-body problem we derived in the previous section 
seems to make no exception from this rule. In what follows, we will point out 
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the singularities that occur in our problem and prove some results related to 
them. The most surprising seems to be the existence of a class of solutions 
with some hybrid singularities, which are both coUisional and non-collisional. 

4.1. Singularities of the equations. The equations of motion flT^ have 
restrictions. First, the variables are constrained to a surface of constant cur- 
vature, i.e. (q, p) G T*(M^)", where M^ is the surface of curvature k 7^ (in 
particular, Mj = S^ and M^^^ = H^), T*(M2)" is the cotangent bundle of 
M^. Second, system (fill) , which contains the gradient ([8]), is undefined in the 
set A := Ui<i<j<.„Aij, with 

A,, :={qG(M2,r |(Kq,0q,.)2 = l}, 

where both the force function ([1]) and its gradient ([8]) become infinite. Thus 
the set A contains the singularities of the equations of motion. 

The singularity condition, [nq^t cy)'^ = 1, suggests that we consider two 
cases, and thus write Aj^ = A^ U A^^-, where 

A+ := {q e (M^)" I «:q,0q^. = 1} and Ar. := {q e (M^)" | «:q,0q,. = -1}. 

Accordingly, we define 

A+ := Ui<i<^<„A,:^ and A" := Ui<,<i<nA-.. 

Then A = A''" U A~. The elements of A+ correspond to collisions for any 
K 7^ 0, whereas the elements of A~ correspond to what we will call antipodal 
singularities when k > 0. The latter occur when two bodies are at the opposite 
ends of the same diameter of a sphere. For k < 0, such singularities do not 
exist because Kqj q^ > 1 . 

In conclusion, the equations of motion are undefined for configurations that 
involve collisions on spheres or hyperboloids, as well as for configurations with 
antipodal bodies on spheres of any curvature k > 0. In both cases, the gravi- 
tational forces become infinite. 

In the 2-body problem. A"*" and A~ are disjoint sets. Indeed, since there 
are only two bodies, /tqi ■ q2 is either +1 or —1, but cannot be both. The 
set A"*" n A~, however, is not empty for n > 3. In the 3-body problem, for 
instance, the configuration in which two bodies are at collision and the third 
lies at the opposite end of the corresponding diameter is, what we will call 
from now on, a collision-antipodal singularity. 

The theory of differential equations merely regards singularities as points 
for which the equations break down, and must therefore be avoided. But 
singularities exhibit sometimes a dynamical structure. In the 3-body problem 
in R, for instance, the set of binary collisions is attractive in the sense that for 
any given initial velocities, there are initial positions such that if two bodies 
come close enough to each other but far enough from other collisions, then the 
collision will take place. (Things are more complicated with triple collisions. 
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Two of the bodies coming close to triple collisions may form a binary while 
the third gets expelled with high velocity away from the other two, [5B].) 

Something similar happens for binary collisions in the 3-body problem on a 
geodesic of S^. Given some initial velocities, one can choose initial positions 
that put mi and m2 close enough to a binary collision, and rris far enough from 
an antipodal singularity with either mi or m,2, such that the binary collision 
takes place. This is indeed the case because the attraction between mi and 
7712 can be made as large as desired by placing the bodies close enough to each 
other. Since m,^ is far enough from an antipodal position, and no comparable 
force can oppose the attraction between mi and m,2, these bodies will collide. 

Antipodal singularities lead to a new phenomenon on geodesies of S^. Given 
initial velocities, no matter how close one chooses initial positions near an 
antipodal singularity, the corresponding solution is repelled in future time 
from this singularity as long as no collision force compensates for this force. 
So while binary collisions can be regarded as attractive if far away from binary 
antipodal singularities, binary antipodal singularities can be seen as repulsive 
if far away from collisions. But what happens when collision and antipodal 
singularities are close to each other? As we will see in the next subsection, 
the behavior of solutions in that region is sensitive to the choice of masses and 
initial conditions. In particular, we will prove the existence of some hybrid 
singular solutions in the 3-body problem, namely those that end in finite time 
in a collision-antipodal singularity, as well as of solutions that reach a collision- 
antipodal configuration but remain analytic at this point. 

4.2. Solution singularities. The set A is related to singularities which arise 
from the question of existence and uniqueness of initial value problems. For 
initial conditions (q, p)(0) G T*(M^)" with q(0) ^ A, standard results of the 
theory of differential equations ensure local existence and uniqueness of an 
analytic solution (q, p) defined on some interval [0, t"*"). Since the surfaces M^ 
are connected, this solution can be analytically extended to an interval [0,t*), 
with < t"*" < t* < oo. If t* = oo, the solution is globally defined. But if 
t* < oo, the solution is called singular, and we say that it has a singularity at 
time t*. 

There is a close connection between singular solutions and singularities of 
the equations of motion. In the classical case {k = 0), this connection was 
pointed out by Paul Painleve towards the end of the 19th century. In his fa- 
mous lectures given in Stockholm, [57] , he showed that every singular solution 
(q, p) is such that q(t) -^ A when t —* t*, for otherwise the solution would 
be globally defined. In the Euclidean case, k = 0, the set A is formed by all 
configurations with collisions, so when q(t) tends to an element of A, the solu- 
tion ends in a collision singularity. But it is also possible that q(t) tends to A 
without asymptotic phase, i.e. by oscillating among various elements without 
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Figure 1 . The relative positions of the force acting on m, while 
the body is on the geodesic z = 0. 

ever reaching a definite position. Painleve conjectured that such noncoUision 
singularities, which he called pseudocoUisions, exist. In 1908, Hugo von Zeipel 
showed that a necessary condition for a solution to experience a pseudocolli- 
sion is that the motion becomes unbounded in finite time, [7S], [SI]. Zhihong 
(Jeff) Xia produced the first example of this kind in 1992, [77j. Historical 
accounts of this development appear in [TH] and [22] . 

The results of Painleve don't remain intact in our problem, [23], [21], so 
whether pseudocoUisions exist for k ^ is not clear. Nevertheless, we will 
now show that there are solutions ending in collision-antipodal singularities of 
the equations of motion, solutions these singularities repel, as well as solutions 
that are not singular at such configurations. To prove these facts, we need the 
result stated below, which provides a criterion for determining the direction of 
motion along a great circle in the framework of an isosceles problem defined 
in an invariant set S^. 



and assume that a body of 
within its first quadrant, 



< 0, the force pulls the body along the circle toward 
> 0, the force pulls the body along the circle toward 



Lemma 1. Consider the n-body problem in S^ 
mass m is at rest at time to on the geodesic z - 
x,y > 0. Then, if 

(a) x(to) > and y(tQ 
the point (x, y) = (1,0). 

(b) x{to) < and y{to 
the point (x, y) = (0, 1) . 

(c) x{to) < and y{to) < 0, the force pulls the body toward the point (1, 0) 
ifyito)/x{to) > y{to)/x{to), toward (0,1) if y{to)/x{tQ) < y{to)/x{to), but no 
force acts on the body if neither of the previous inequalities holds. 

(d) x(to) > and y{tQ) > 0, the motion is impossible. 

Proof. By equation (Q, xx + yy = — (x^ + y"^) < 0, which means that the 
force acting on m is always directed along the tangent at m to the geodesic 
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circle 2; = or inside the half-plane containing this circle. Assuming that an 
xy-coordinate system is fixed at the origin of the acceleration vector (point P 
in Figure [1]), this vector always lies in the half-plane below the line of slope 
—xitoj/yito) (i.e. the tangent to the circle at the point P in Figure [1]). We 
further prove each case separately. 

(a) If xito) > and i/(to) < 0, the force acting on m is represented by 
a vector that lies in the region given by the intersection of the fourth quad- 
rant (counted counterclockwise) and the half plane below the line of slope 
— x(to)/|/(^o)- Then, obviously, the force pulls the body along the circle in the 
direction of the point (1, 0). 

(b) If x(to) < and y(to) > 0, the force acting on m is represented by a 
vector that lies in the region given by the intersection of the second quadrant 
and the half plane lying below the line of slope — x(to)/y(^o)- Then, obviously, 
the force pulls the body along the circle in the direction of the point (0, 1). 

(c) If x(to) < and y(to) < 0, the force acting on m is represented by 
a vector lying in the third quadrant. Then the direction in which this force 
acts depends on whether the acceleration vector lies: (i) below the line of 
slope y(to)/x(to) (PB is below OP in Figure [T]); (ii) above it (PC is above 
OP); or (iii) on it (i.e. on the line OP). Case (iii) includes the case when the 
acceleration is zero. 

In case (i), the acceleration vector lies on a line whose slope is larger than 
y(tQ)/x(to), i.e. y(to)/x(tQ) > y(to)/x(tQ), so the force pulls m toward (1,0). 
In case (ii) , the acceleration vector lies on a line of slope that is smaller than 
y{to)/x{to), i.e. y(to)/x(tQ) < y(to)/x{to), so the force pulls m toward (0,1). 
In case (iii), the acceleration vector is either zero or lies on the line of slope 
y(tQ)/x(to), i.e. i/(to)/x(to) = y{tQ)/x(to). But the latter alternative never 
happens. This fact follows from the equations of motion ([7]), which show that 
the acceleration is the difference between the gradient of the force function 
and a multiple of the position vector. But according to Euler's formula for 
homogeneous functions, ([3]), and the fact that the velocities are zero, these 
vectors are orthogonal, so their difference can have the same direction as one 
of them only if it is zero. This vectorial argument agrees with the kinematic 
facts, which show that if x(to) = y{to) = and the acceleration has the same 
direction as the position vector, then m doesn't move, so x{t) = y{t) = 0, 
and therefore x{t) = y{t) = for all t. In particular, this means that when 
i/(to) = x(to) = 0, no force acts on m, so the body remains fixed. 

(d) If x(to) > and y(to) > 0, the force acting on m is represented by 
a vector that lies in the region given by the intersection between the first 
quadrant and the half-plane lying below the line of slope —a; (to)/?/ (to). But 
this region is empty, so the motion doesn't take place. D 
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m^ =: m 
Figure 2. The initial positions of mi, 

desic z = 0. 



7712, and 777,3 on the geo- 



We will further prove the existence of solutions with collision-antipodal sin- 
gularities, solutions repelled from collision-antipodal singularities in positive 
time, as well as of solutions that remain analytic at a collision-antipodal con- 
figuration. They show that the dynamics of A"*" fl A~ is more complicated 
than the dynamics of A"*" and A~ away from the intersection, since solutions 
can go both towards and away from this set for t > 0, and can even avoid sin- 
gularities. This result represents a first example of a non-collision singularity 
reached by only three bodies as well as a first example of a non-singularity 
collision. 

Theorem 1. Co7isider the 3-body proble7n m S^ with the bodies mi a7id ?7i2 
havi7ig 7nass M > aTid the body m^ haviTig 7nass m > 0. TheTi 

(i) there are values of m and M, as well as initial conditions, for which the 
solutions end in finite time in a collision- antipodal singularity; 

(a) other choices of masses and initial conditions lead to solutions that are 
repelled from a collision- antipodal singularity; 

(Hi) and yet other choices of masses and initial data correspond to solutions 
that reach a collision- antipodal configuration but remain analytic at this point. 

Proof. Let us start with some initial conditions we will refine on the way. 
During the refinement process, we will also choose suitable masses. Consider 

xi(0) = -x(0), yi(0) = 2/(0), ^i(0)=0, 

X2(0) = a;(0), y2{0) = y{0), ^2(0) = 0, 

X3(0) = 0, ^3(0) = -1, ^3(0) = 0, 

as well as zero initial velocities, where < x{t),y(t) < 1 are functions with 
x{ty + |/(t)^ = 1. Since all z coordinates are zero, only the equations of 
coordinates x and y play a role in the motion. The symmetry of these initial 
conditions implies that m-s remains fixed for all time (in fact the equations 
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corresponding to X3 and jjs reduce to identities), that the angular momentum 
is zero, and that it is enough to see what happens for m2, because mi behaves 
symmetrically with respect to the y axis. Thus, substituting the above initial 
conditions into the equations of motion, we obtain 

These equations show that several situations occur, depending on the choice 
of masses and initial positions. Here are two significant possibilities. 

1. For M > Am, it follows that x(0) < and y{0) > for any choices of 
initial positions with < x(0), y{0) < 1. 

2. For M < 4m, there are initial positions for which: 

(a) x(0) < and y{0) > 0, 

(b) x(0) > and i/(0) < 0, 

(c) x(0) = m = 0. 

In case 2(c), the solutions are fixed points of the equations of motion, a 
situation achieved, for instance, when M = 2m, and a;(0) = y(0) = \/2/2. The 
cases of interest for us, however, are 1 and 2(b). In the former, 7712 begins to 
move from rest towards a collision with mi at (0, 1), but whether this collision 
takes place also depends on velocities, which affect the equations of motion. 
In the latter case, m2 moves away from the same collision, and we need to 
see again how the velocities alter this initial tendency. So let us write now 
the equations of motion for m2 starting from arbitrary masses M and m. The 
computations lead us to the system 

^ ^ \y = ^-^-ii' + y')y 

and the energy integral 

.2^.2 h 2my M{2y^-l) 
X +y = -— h 



M X 2xy 

Substituting this expression of x"^ + y"^ into equations (!34l) . we obtain 

{- _ 4{M-2m)x'^-2{M-2m)x^-M+4m h_ 
^ ~ 4x2?; M^ 

■■ ^ M+2{M-2m)y^-A{M-2m)y'i _ h 
y 4xj/2 My- 

We will further focus on the first class of orbits announced in this theorem. 

(i) To prove the existence of solutions with collision-antipodal singularities, 
let us further examine the case M = 8m, which brings system (l35l) to the form 

{" 6mx^ 3m m h 
y y x^y 8m 

y = 2m _|_ 3m _ 6my^ h_ 
" xy"^ X X 8m"' 
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with the energy integral 

Amx 2my h 



(37) x' + y' + 



y X am 



Then, as x — ;► and ?/ — ;► 1, both x and ij tend to — oo, so they are ulti- 
mately negative. This fact corresponds to case (c) of Lemma [H But a simple 
computation shows that y/x tends to zero as x — *> and y — > 1. Since 
y/x > 0, it follows that if {x{0),y{0)) is chosen close enough to (0,1), then 
i/(0)/x(0) < y{0)/x{0), so according to the conclusion of Lemma [T](c) the 
collision-antipodal configuration is reached. As the forces and the potential 
are infinite at this point, using the energy relation (1371) it follows that the 
velocities are also infinite. Consequently the motion cannot be analytically 
extended beyond the collision-antipodal configuration, which thus proves to 
be a singularity. 

(ii) To show the existence of solutions repelled from a collision-antipodal 
singularity of the equations of motion in positive time, let us take M = 2m. 
Then equations fl35|) have the form 

{^ — m h_^ 
y 2xy'2 2m y' 

with the integral of energy 

(39) x' + f + ^ = ^, 

xy 2m 

which implies that h > 0. Obviously, as x — > and y —^ 1, the forces and and 
the kinetic energy become infinite, so the collision-antipodal configuration is 
a singularity if it were to be reached. But as we will further see, this can- 
not happen for this choice of masses. Indeed, as we saw in case 2(c) above, 
the initial position x(0) = y{0) = \/2/2 corresponds to a fixed point of the 
equations of motion for zero initial velocities. Therefore we must seek the 
desired solution for initial conditions with < x(0) < V^/2 and the corre- 
sponding choice of 1/(0) > 0. Let us pick any such initial positions, as close to 
the collision-antipodal singularity as we want, and zero initial velocities. For 
X — > 0, however, equations ( l38l) show that both x and y grow positive. But 
according to case (d) of Lemma [H such an outcome is impossible, so the mo- 
tion cannot come infinitesimally close to the corresponding collision-antipodal 
singularity, which repels any solution with M = 2m and initial conditions 
chosen as we previously described. 

(iii) To prove the existence of solutions that have no singularity at a collision- 
antipodal configuration, let us further examine the case M = 4m, which brings 
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system ( 135|) to the form 



(40) 



y Am 

mx(2y^+l) h 

y"^ 4m ^ ' 



For this choice of masses, the energy integral becomes 

,^_,. •2-2 2mx h 

(41) x^ + y^ + = — . 

y Am 

We can compute the value of h from the initial conditions. Thus, for initial 
positions x(0),y(0) and initial velocities x(0) = ?/(0) = 0, the energy constant 
is /i = 8m2x(0)/y(0) > 0. 

Assuming that x — ^ and y — i> 1, equations (l40l) imply that x{t) —* —m < 
and y{t) -^ —h/Am < 0, which means that the forces are finite at the 
collision-antipodal configuration. We are thus in the case (c) of Lemma [H so 
to determine the direction of motion for 7712 when it comes close to (0, 1), we 
need to take into account the ratio y/x, which tends to h/Am? as x — > 0. Since 
h = 8m^x{0)/y{0), lim^_o(i//x) = 2x{0)/y{0). Then 2x{0)/y{0) < y{0)/x{0) 
for any x(0) and y{0) with < x(0) < l/VS and the corresponding choice 
of 1/(0) > given by the constraint x^(0) + y^(0) = 1. But the inequality 
2x{0)/y{0) < y{0)/x{0) is equivalent to the condition i/(to)/^(^o) < yito)/x(to) 
in Lemma [T](c), according to which the force pulls 777,2 toward (0, 1). Therefore 
the velocity and the force acting on m2 keep this body on the same path until 
the collision- antipodal configuration occurs. 

It is also clear from equation fHT]) that the velocity is positive and finite 
at collision. Since the distance between the initial position and (0, 1) is 
bounded, 7712 collides with mi in finite time. Therefore the choice of masses 
with M = 4??7,, initial positions x{0),y{0) with < x(0) < l/v3 and the 
corresponding value of y{0), and initial velocities x{0) = y{0) = 0, leads to a 
solution that remains analytic at the collision-antipodal configuration, so the 
motion naturally extends beyond this point. D 

5. Relative equilibria in S^ 

In this section we will prove a few results related to fixed points and elliptic 
relative equilibria in S^. Since, by Euler's theorem (see Appendix), every 
element of the group 5*0(3) can be written, in an orthonormal basis, as a 
rotation about the z axis, we can define elliptic relative equilibria as follows. 

Definition 1. An elliptic relative equilibrium in S^ is a solution of the form 
Qj = {xi, yi,Zi), i = 1, . . . , n, of equations /[W\) with Xi = ri cos(ci;t + ai),yi = 
Ti sin(ci;t -|- ctj), Zi = constant, where u, ai, and Vi, with < rj = (1 — zfY^'^ < 
1, i = 1, . . . ,n, are constants. 
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Notice that although the equations of motion don't have an integral of 
the center of mass, a "weak" property of this kind occurs for elliptic relative 
equilibria. Indeed, it is easy to see that if all the bodies are at all times on one 
side of a plane containing the rotation axis, then the integrals of the angular 
momentum are violated. This happens because under such circumstances the 
vector representing the total angular momentum cannot be zero or parallel to 
the z axis. 

5.1. Fixed points. The simplest solutions of the equations of motion are 
fixed points. They can be seen as trivial relative equilibria that correspond to 
cj = 0. In terms of the equations of motion, we can define them as follows. 

Definition 2. A solution of system ^2^) is called a fixed point if 

Vq, t/i(q)(t) = pi(t) = for all t G R and i = 1, . . . , n. 

Let us start with finding the simplest fixed points, namely those that occur 
when all the masses are equal. 

Theorem 2. Consider the n-body problem in S^ with n odd. If the masses 
are all equal, the regular n-gon lying on any geodesic is a fixed point of the 
equations of motion. For n = 4, the regular tetrahedron is a fixed point too. 

Proof. Assume that nii = m2 = ■ ■ ■ = nin, and consider an n-gon with an odd 
number of sides inscribed in a geodesic of S^ with a body, initially at rest, 
at each vertex. In general, two forces act on the body of mass mf. the force 
Vqi?7i(q), which is due to the interaction with the other bodies, and the force 
—mi{c[i ■ qi)qi, which is due to the constraints. The latter force is zero at t = 
because the bodies are initially at rest. Since q^ ■ Vq.f/i(q) = 0, it follows that 
Vq^f/i(q) is orthogonal to q^, and thus tangent to S^. Then the symmetry of 
the n-gon implies that, at the initial moment t = 0, Vq;f/i(q) is the sum of 
pairs of forces, each pair consisting of opposite forces that cancel each other. 
This means that Vq.f/i(q) = 0. Therefore, from the equations of motion and 
the fact that the bodies are initially at rest, it follows that 

qi(0) = -(q.(0) ■ q.(0))qi(0) = 0, z = 1, . . . , n. 

But then no force acts on the body of mass nii at time t = 0, consequently 
the body doesn't move. So qi(t) = for all t G R. Then qi(t) = for all 
t G R, therefore Vq.f/i(q)(t) = for all t G R, so the n-gon is a fixed point 
of equations (ITI?]) . 

Notice that if n is even, the n-gon has n/2 pairs of antipodal vertices. Since 
antipodal bodies introduce singularities into the equations of motion, only the 
n-gons with an odd number of vertices are fixed points of equations (TT9|) . 

The proof that the regular tetrahedron is a fixed point can be merely done 
by computing that 4 bodies of equal masses with initial coordinates given 
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by qi = (0,0,l),q2 = (0, 2v^/3, -1/3), qa = (-2/v^, -v^/3, -1/3), q4 = 
(2/\/6, — \/2/3, —1/3), satisfy system flT9l) . or by noticing that the forces acting 
on each body cancel each other because of the involved symmetry. D 

Remark 1. If equal masses are placed at the vertices of the other four regular 
polyhedra: octahedron (6 bodies), cube (8 bodies), dodecahedron (12 bodies), 
and icosahedron (20 bodies), they do not form fixed points because antipodal 
singularities occur in each case. 

Remark 2. In the proof of Theorem [H we discovered that if one body has 
mass m and the other two mass M = 2m, then the isosceles triangle with 
the vertices at (0,-1,0), (-^2/2, \/2/2, 0), and (^2/2,^2/2,0) is a fixed 
point. Therefore one might expect that fixed points can be found for any 
given masses. But, as formula ( !33|) shows, this is not the case. Indeed, if 
one body has mass m and the other two have masses M > Am, there is no 
configuration (which must be isosceles due to symmetry) that corresponds to 
a fixed point since x and y are never zero. This observation proves that in the 
3-body problem, there are choices of masses for which the equations of motion 
lack fixed points. 

The following statement is an obvious consequence of the proof given for 
Theorem [2l 

Corollary 1. Consider an odd number of equal bodies, initially at the ver- 
tices of a regular n-gon inscribed in a great circle of S"^ , and assume that the 
solution generated from this initial position maintains the same relative config- 
uration for all times. Then, for all t G R, this solution satisfies the conditions 
Vq,f/i(q(t)) = 0, i = l,...,n. 

It is interesting to see that if the bodies are within a hemisphere (meaning 
half a sphere and its geodesic boundary), fixed points do not occur if at least 
one body is not on the boundary. Let us formally state and prove this result. 

Theorem 3. Consider an initial nonsingular configuration of the n-body prob- 
lem in S^ for which all bodies lie within a hemisphere, meant to include its 
geodesic boundary, with at least one body not on this geodesic. Then this con- 
figuration is not a fixed point. 

Proof. Without loss of generality we can consider the initial configuration of 
the bodies m,i, . . . , rrin in the hemisphere z > 0, whose boundary is the geodesic 
z = 0. Then at least one body has the smallest z coordinate, and let mi be one 
of these bodies. Also, at least one body has its z coordinate positive, and let 
7712 be one of them. Since all initial velocities are zero, only the mutual forces 
between bodies act on mi. Then, according to the equations of motion (fTTll . 
^^i-2i(0) = ^f/i(q(0)). But as no body has its z coordinate smaller than Zi, 



The n-Body Problem in Spaces of Constant Curvature 33 

the terms contained in the expression of ^t/i(q(0)) that involve interactions 
between mi and rrii are all larger than or equal to zero for z = 3, 4, . . . , ra, while 
the term involving m2 is strictly positive. Therefore ^f/i(q(0)) > 0, so nii 
moves upward the hemisphere. Consequently the initial configuration is not a 
fixed point. D 

5.2. Polygonal solutions. We will further show that fixed points lying on 
geodesies of spheres can generate relative equilibria. 

Theorem 4. Consider a fixed point given by the masses mi, m,2, . . . , m,n that 
lie on a great circle of S^ . Then for every nonzero angular velocity, this con- 
figuration generates a relative equilibrium along the great circle. 

Proof. Without loss of generality, we assume that the great circle is the equator 
z = and that for some given masses mi, m2, . . . , m„ there exist ai, 0^2, . . . , an 
such that the configuration q = (qi,...,qn) given by q^ = {xi,yi,0),i = 
1, . . . ,n, with 

(42) Xi = cos{ujt + ttj), yi = sin(u;t + aj), i = 1, . . . ,n, 

is a fixed point for u = 0. This configuration can also be interpreted as being 
q(0), i.e. the solution q at t = for any cu 7^ 0. So we can conclude that 
Vq-f/i(q(0)) = 0, i = 1, . . . ,n. But then, for t = 0, the equations of motion 



Xi = -{xj + yf)xi 



fITTj) reduce to 
(43) 



i = 1, . . . ,n. Notice, however, that Xi = — ct;sin(a;t + ai),Xi = —lj^ cos{ujt + 
ttj)) Vi = —^cos{ujt + ai), and jji = —uP' sin(a;t + a^), therefore x\ + ijl = uj^ . 
Using these computations, it is easy to see that q given by (1121) is a solution of 
(H5]l for every t, so no forces due to the constraints act on the bodies, neither 
at t = nor later. Since Vq.f/i(q(0)) = 0, z = 1, . . . ,?7,, it follows that the 
gravitational forces are in equilibrium at the initial moment, so no gravitational 
forces act on the bodies either. Consequently, the rotation imposed by a; 7^ 
makes the system move like a rigid body, so the gravitational forces further 
remain in equilibrium, consequently Vq,t/i(q(t)) = 0, i = 1, . . . ,n, for all t. 
Therefore q given by (H2l) satisfies equations (TT7|) . Then, by Definition [1], q is 
an elliptic relative equilibrium. D 

The following result shows that relative equilibria generated by fixed points 
obtained from regular n-gons on a great circle of S^ can occur only when the 
bodies rotate along the great circle. 

Theorem 5. Consider an odd number of equal bodies, initially at the vertices 
of a regular n-gon inscribed in a great circle of S^. Then the only elliptic 
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relative equilibria that can be generated from this configuration are the ones 
that rotate in the plane of the original great circle. 

Proof. Without loss of generality, we can prove this result for the equator 
z = 0. Consider therefore an elliptic relative equilibrium solution of the form 

(44) Xi = Ti cos{ujt + Oj), Hi = Ti sm{u)t + «»), Zi = ±{l ~ r,^)^/^ 

i = 1, . . . , n, with + taken for 2, > and — for z.i < 0. The only condition we 
impose on this solution is that r^ and ctj, i = 1, . . . ,n, are chosen such that 
the configuration is a regular n-gon inscribed in a moving great circle of S^ at 
all times. Therefore the plane of the n-gon can have any angle with, say, the 
z-axis. This solution has the derivatives 

Xi = —riUJ?,m.{ujt + ai), yi = riU cos{ujt + ai), ij = 0, i = 1, . . . ,n, 
Xi = —TiUJ^ cos(ci;t + OLi), xji = —TiU? sin(u;t + ttj), Zi = 0, i = 1, . . . ,n. 



Then 



Xi+yi+Zi=riUJ,i = l,...,n. 



Since, by Corollary [H any n-gon solution with n odd satisfies the conditions 

Vq,f/i(q)=0, z = l,...,n, 
system flT9|) reduces to 



= -{il + yf + zf)xi, 

-- -(if + yf + zf)yi, 

-- -{il + yf + zf)zi, i = l,...,n. 
Then the substitution of (jUj) into the above equations leads to: 

rj(l — rf)uj^ cos{u!t + a,) = 0, 

rj(l — r^)iu'^ sin(u;t + Oj) = 0, i = 1, . . . , n. 

But assuming u ^ 0, this system is nontrivially satisfied if and only if Tj = 1, 
conditions which are equivalent to 2;j = 0, i = 1, . . . ,n. Therefore the bodies 
must rotate along the equator z = 0. D 

Theorem [5] raises the question whether elliptic relative equilibria given by 
regular polygons can rotate on other curves than geodesies. The answer is 
given by the following result. 

Theorem 6. Consider the n-body problem with equal masses in S^. Then, for 
any n odd, m > and z G (—1, 1), there are a positive and a negative u that 
produce elliptic relative equilibria in which the bodies are at the vertices of an 
n-gon rotating in the plane z = constant. If n is even, this property is still 
true if we exclude the case z = 0. 



The n-Body Problem in Spaces of Constant Curvature 35 

Proof. There are two cases to discuss: (i) n odd and (ii) n even. 

(i) To simplify the presentation, we further denote the bodies by mi,i = 
—s, — s + 1, . . . , —1, 0, 1, . . . , s — 1, s, where s is a positive integer, and as- 
sume that they all have mass m. Without loss of generality we can fur- 
ther substitute into equations (fT9l) a solution of the form (jUj) with i as 

1 4s7r Stt r\ Stt 4s7r 

auOVe, Cl—s Is+T' ■ ■ ■ ' ^—l 2s+l ' '^O ' ^^ 2s+l ' ■ ■ ■ ' ^^ 2s+l ' 

r := ri,z := Zi, and consider only the equations for i = 0. The study of this 
case suffices due to the involved symmetry, which yields the same conclusions 
for any value of i. 

The equation corresponding to the zq coordinate takes the form 

Y^ m^z — kojz) 2 2 



-±1.n (1 - kl)'/' 



r^u^z = 0, 



where koj = xqXj + yoyj + zqZj = cosaj — z"^ cosa^ + z"^. Using the fact that 
Yj = cosa_j, and /cqj = kQ^^-j), 

,2 



j,2 _|_ ^2 _ -^^ cos a,- = cosa_o, and k^j = kQ/j), this equation becomes 



E^ 



cosajj uj 



Now we need to check whether the equations corresponding to Xq and i/q lead 
to the same equation. In fact, checking for xq, and ignoring i/q, suffices due to 
the same symmetry reasons invoked earlier or the duality of the trigonometric 
functions sin and cos. The substitution of the the above functions into the 
first equation of (fT9!) leads us to 

2 ^^ 2 J. \;^ m[cos{ujt + aj) — koj cosLut] 



l)uj'^ coscut = y. 



A straightforward computation, which uses the fact that r"^ + z"^ = 1, sina^ = 
— sina_j, costtj = cosa_j, and k^j = ko(^_j), yields the same equation fHSl) . 
Writing the denominator of equation (H5]l explicitly, we are led to 



(46) 



^ (i - cos a,- )"/"(! - z'']-''''\2 - (i - COS a,) (i - z^]]"'" m 



-^{1- cosa,)V2(i _ ^2)3/2[2 _ (1 _ cosa,)(l - z^)f/^ 



The left hand side is always positive, so for any ?7i > and z G (—1,1) fixed, 
there are a positive and a negative u that satisfy the equation. Therefore the 
n-gon with an odd number of sides is an elliptic relative equilibrium. 

(ii) To simplify the presentation when n is even, we denote the bodies 
by rrii, i = —s + 1,..., —1, 0, 1, . . . , s — 1, s, where s is a positive integer, 
and assume that they all have mass m. Without loss of generality, we can 
substitute into equations 0191) a solution of the form 044p with i as above, 

( — •5 + 1)71" TT n TT (s— l)7r 

a_s+i = - — r^, . . . , a_i = --, ao = 0, ai = -,..., a^-i = ^-7^, "s = vr. 
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r := Ti, z := Zi, and consider as in the previous case only the equations for 
i = 0. Then using the fact that k^j = /co(-j), cosa^ = cosa_j, and cosvr = — 1, 
a straightforward computation brings the equation corresponding to Zq to the 
form 



(47) y^ 2(1 -cos a,- 



kiy/' {i-kiy/' 



co^ 



ra 



Using additionally the relations sinoj = — sina.j and sinvr = 0, we obtain 
for the equation corresponding to xq the same form (1471) . which — for koj and 
kos written explicitly — becomes 



s-l 



^n 



.^^ ,_ cosaj)V2(i _ ^2)3/2[2 _ (1 _ cosaj){l - z^)]^^ 

+ 



iU^ 



4z^\z\{l - z^f/^ m' 

Since the left hand side of this equations is positive and finite, given any ?Ti > 
and z G (—1, 0) U (0, 1), there are a positive and a negative a; that satisfy it. So 
except for the case z = 0, which introduces antipodal singularities, the rotating 
ra-gon with an even number of sides is an elliptic relative equilibrium. D 

5.3. Lagrangian solutions. The case n = 3 presents particular interest in 
the Euclidean case because the equilateral triangle is an elliptic relative equi- 
librium for any values of the masses, not only when the masses are equal. But 
before we check whether this fact holds in S^, let us consider the case of three 
equal masses in more detail. 

Corollary 2. Consider the 3-body problem with equal masses, m := mi = 
7722 = rris, in S^. Then for any m > and z G (—1, 1), there are a positive 
and a negative uo that produce elliptic relative equilibria in which the bodies 
are at the vertices of an equilateral triangle that rotates in the plane z = 
constant. Moreover, for every u'^/m there are two values of z that lead to 
relative equilibria if uj'^/m E (8/-\/3, oo) U {3}, three values if uj'^/m = 8/a/3, 
and four values if ij^jm E (3,8/-\/3). 

Proof. The first part of the statement is a consequence of TheoremElfor n = 3. 
Alternatively, we can substitute into system (IT9l) a solution of the form (jHl) 
with 2 = 1, 2, 3, r := ri = r2 = r^, z = ±(1 — r^)^/^, ai = 0, a2 = 27r/3, as = 
47r/3, and obtain the equation 

(48) ^ ? -^^. 

^ ^ V3(l + 2^2 _ 3^4)3/2 rn 
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Figure 3. The graph of the function f{z) = ^ ^ ^ 3^^ for 
ZG(-1,1). 



The left hand side is positive for z e (—1,1) and tends to infinity when z —>■ ±1 
(see Figure [3l). So for any z in this interval and m > 0, there are a positive 
and a negative u for which the above equation is satisfied. Figure [3] and a 
straightforward computation also clarify the second part of the statement. D 



Remark 3. A result similar to Corollary [2] can be proved for two equal masses 
that rotate on a non-geodesic circle, when the bodies are situated at opposite 
ends of a rotating diameter. Then, for z G (—1,0) U (0, 1), the analogue of 
( l48l) is the equation 



ui 



m 



The case z 



4z2|^|(l-^2)3/2 

yields no solution because it involves an antipodal singularity. 



We have reached now the point when we can decide whether the equilateral 
triangle can be an elliptic relative equilibrium in S^ if the masses are not equal. 
The following result shows that, unlike in the Euclidean case, the answer is 
negative when the bodies move on the sphere in the same Euclidean plane. 

Proposition 2. In the 3-body problem in S^, if the bodies 7711,7712,1^3 are 
initially at the vertices of an equilateral triangle in the plane z = constant 
for some z G (—1,1), then there are initial velocities that lead to an elliptic 
relative equilibrium in which the triangle rotates in its own plane if and only 

if TTll = 777,2 = "^3- 

Proof. The implication which shows that if mi = m2 = m^, the rotating 
equilateral triangle is a relative equilibrium, follows from Theorem [21 To 
prove the other implication, we substitute into equations (HM a solution of the 
form (144]) with i = 1,2,3, r := ri,r2,rs, z := Zi = Z2 = z^ = ±(1 — r'^Y^'^, 
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and «! = 0, 02 = 27r/3, as = 47r/3. The computations then lead to the system 

{mi + m2 = '^uP' 
1712 + ^^3 = 1^"^ 
ms + mi = 7u;^, 

where 7 = a/3(1 + 22;^ — 32^)^/^/4. But for any z = constant in the interval 
(—1, 1), the above system has a solution only for nii = m2 = m^ = 'yuj'^/2. 
Therefore the masses must be equal. D 

The next result leads to the conclusion that Lagrangian solutions in S^ can 
take place only in Euclidean planes of R^. This property is known to be true 
in the Euclidean case for all elliptic relative equilibria, [76], but Wintner's 
proof doesn't work in our case because it uses the integral of the center of 
mass. Most importantly, our result also implies that Lagrangian orbits with 
non-equal masses cannot exist in S^. 

Theorem 7. For all Lagrangian solutions in S^, the masses mi,m2 and m^ 
have to rotate on the same circle, whose plane must he orthogonal to the rota- 
tion axis, and therefore mi = m2 = m-^ . 

Proof. Consider a Lagrangian solution in S^ with bodies of masses mi,m2, 
and 7713. This means that the solution, which is an elliptic relative equilibrium, 
must have the form 



Xi = Ti cos cut, yi = Ti sin cut, 2:1 = (1 — r^ 

X2 = r2 cos{ujt + a), 2/2 = ^"2 sin(u;t + a), Z2 = {I — r. 



2N1/2 



) 



iy^\ 



X3 = rg cos{ujt + b), y3 = r-i sm{u;t + b), 2:3 = (1 - rlY^^, 

with 6 > a > 0. In other words, we assume that this equilateral forms a 
constant angle with the rotation axis, z, such that each body describes its 
own circle on S^. But for such a solution to exist, it is necessary that the total 
angular momentum is either zero or is given by a vector parallel with the z axis. 
Otherwise this vector rotates around the z axis, in violation of the angular- 
momentum integrals. This means that at least the first two components of 
the vector Yli=i "^jQi ^ Qi must be zero. A straightforward computation shows 
this constraint to lead to the condition 

niiTiZi sin ut + ?7i2?"2^2 sin(a;t + a) + m^r^z^ sin(a;t + b) = 0, 

assuming that u ^ 0. For t = 0, this equation becomes 

(50) 777,2^222 siu a = —m-^r^z^ sin b. 

Using now the fact that 

a := X1X2 + yiy2 + 21^2 = 3:1X3 + yiy^ + 2123 = X3X2 + 1/32/2 + Z3Z2 
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is constant because the triangle is equilateral, the equation of the system of 
motion corresponding to jji takes the form 

Kri(rl — l)uj'^ sincut = m2^2 sm(ujt + a) + rrisrs sin(ci;t + b), 

where K is a. nonzero constant. For t = 0, this equation becomes 

(51) 777,2^2 sin a = —m^r^ sin h. 

Dividing (l50l) by (15T1) . we obtain that Z2 = Z3. Similarly, we can show that 
Zi = Z2 = z^, therefore the motion must take place in the same Euclidian 
plane on a circle orthogonal to the rotation axis. Proposition [2] then implies 

that mi = 7712 = 1^3- n 

5.4. Eulerian solutions. It is now natural to ask whether such elliptic rel- 
ative equilibria exist, since — as Theorem [5] shows — they cannot be generated 
from regular n-gons. The answer in the case n = 3 of equal masses is given by 
the following result. 

Theorem 8. Co7isider the 3-hody probleTn m S^ with equal 7nasses, m := 7ni = 
^2 = iTT'S- Fix the body of Tnass rrii at (0,0, 1) and the bodies of Tnasses m2 
and ms at the opposite ends of a diameter on the circle z = constant. Then, 
for any m > and z e (—0.5, 0) U (0, 1), there are a positive and a negative u 
that produce elliptic relative equilibria. 

Proof. Substituting into the equations of motion (IT9l) a solution of the form 

xi = 0, yi = 0, zi = 1, 

X2 = r cos tut, y2 = r sin ut, Z2 = z, 

X3 = rcos{ujt + vr), y^ = r sin(co't + vr), 23 = z, 

with r > and z constants satisfying r^ + 2^ = 1, leads either to identities or 
to the algebraic equation 

(52) ^z+\z\-' ^ ^ 

W ) 4^2(1-22)3/2 m' 

The function on the left hand side is negative for z e (—1, —0.5), at 2; = —0.5, 
positive for z G (—0.5, 0) U (0, 1), and undefined at 2; = 0. Therefore, for every 
m > and z G (—0.5, 0) U (0, 1), there are a positive and a negative u that lead 
to a geodesic relative equilibrium. For z = —0.5, we recover the equilateral 
fixed point. The sign of uj determines the sense of rotation. D 

Remark 4. For every u'^/m G (0, 64vT5/45), there are three values of z that 
satisfy relation ( l52l) : one in the interval (—0.5, 0) and two in the interval (0, 1) 
(see Figure Hj). 
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Figure 4. The graph of the function f{z) 
intervals (—1,0) and (0,1), respectively. 



4z+|2|-i 

422(1-^2)3/2 



in the 



Remark 5. If in Theorem|H]we take the masses mi 
the analogue of equation ( l52l) is 

4mz + M\z\~'^ _ u"^ 

4^2(1 - 2;2)3/2 ^ ■ 



=: m and 7722 = m^ =: M, 



Then solutions exist for any z E {—yM/m/2^Qi) U (0,1). This means that 
there are no fixed points for M > Am (a fact that agrees with what we learned 
from Remark [2] and the proof of Theorem [1]) , so relative equilibria exist for 
such masses for all z G (—1, 0) U (0, 1). 



6. Relative equilibria in H^ 

In this section we will prove a few results about fixed points, as well as 
elliptic and hyperbolic relative equilibria in H^. We also show that parabolic 
relative equilibria do not exist. Since, by the Principal Axis theorem for the 
Lorentz group, every Lorentzian rotation (see Appendix) can be written, in 
some basis, either as an elliptic rotation about the z axis, or as an hyperbolic 
rotation about the x axis, or as a parabolic rotation about the line x = 0,y = z, 
we can define three kinds of relative equilibria: the elliptic relative equilibria, 
the hyperbolic relative equilibria, and the parabolic relative equilibria. This 
terminology matches the standard terminology of hyperbolic geometry [53] . 

The elliptic relative equilibria are defined as follows. 

Definition 3. An elliptic relative equilibrium in H^ is a solution q, = [xi, i/i, Zi), 
i = 1, . . . ,n, of equations ([l^j with Xi = pi cos(u;t + a^), yi = pi sin(a;t -|- ai), 
and Zi = [pf + lY^"^ , where to, ai, and pi, i = 1, . . . ,n, are constants. 

Remark that, as in S^, a "weak" property of the center of mass occurs in 
H^ for elliptic relative equilibria. Indeed, if all the bodies are at all times on 
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one side of a plane containing the rotation axis, then the integrals of the angu- 
lar momentum are violated because the vector representing the total angular 
momentum cannot be zero or parallel to the z axis. 

Let us now define the hyperbolic relative equilibria. 

Definition 4. A hyperbolic relative equilibrium in H^ is a solution of equations 
(d^j of the form qj = (xj, yi,Zi), i = 1, . . . , n, defined for all t G R, with 

(53) Xi = constant, i/i = pi sinh(ci;t + Oj), and Zi = pi cosh(ct;t + Oj), 
where u, a,, and pj = (1 + xfY^"^ > 1, i = 1, . . . , n, are constants. 

Finally, the parabolic relative equilibria are defined as follows. 

Definition 5. A parabolic relative equilibrium in H^ is a solution of equations 
(d^j of the form q^ = {xi, yi,Zi), i = 1, . . . , n, defined for all t G R, with 

tJyi — tX^ ^i ^ I L'j fy 

(54) yi = ait + bi{l-ty2) + City2 

Zi = ait-he/2 + c,{l+e/2), 

where ai, bi and Ci, i = 1, . . . , n, are constants, and af + bf — of = —1. 

6.1. Fixed Points in H^. The simplest solutions of the equations of motion 
are the fixed points. They can be seen as trivial elliptic relative equilibria that 
correspond to cu = 0. In terms of the equations of motion, we can define them 
as follows. 

Definition 6. A solution of system [29^) is called a fixed point if 

Vq,f/_i(q)(t) = Pi(t) = for all tGR and z = l,...,ra. 

Unlike in S^, there are no fixed points in H^. Let us formally state and 
prove this fact. 

Theorem 9. In the n-body problem with n > 2 in H^ there are no configura- 
tions that correspond to fixed points of the equations of motion. 

Proof. Consider any coUisionless configuration of n bodies initially at rest in 
H^. This means that the component of the forces acting on bodies due to the 
constraints, which involve the factors xj + yf — zf, i = 1, . . . ,n, are zero at 
t = 0. At least one body, rui, has the largest z coordinate. Notice that the 
interaction between rrii and any other body takes place along geodesies, which 
are concave-up hyperbolas on the {z > 0)-sheet of the hyperboloid modeling 
H^. Then the body rrij^j ^ i, exercises an attraction on rrii down the geodesic 
hyperbola that connects these bodies, so the z coordinate of this force acting 
on nii is negative, independently of whether Zj{0) < Zi{0) or Zj{0) = Zi{0). 
Since this is true for every j = l,...,n, j 7^ i, it follows that Zi{0) < 0. 
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Therefore rrii moves downwards the hyperboloid, so the original configuration 
is not a fixed point. D 

6.2. Elliptic Relative Equilibria in H^. We now consider elliptic relative 
equilibria, and prove an analogue of Theorem [61 

Theorem 10. Consider the n-body problem with equal masses in H^. Then, 
for any m > and z > 1, there are a positive and a negative u that produce 
elliptic relative equilibria in which the bodies are at the vertices of an n-gon 
rotating in the plane z = constant. 

Proof. The proof works in the same way as for Theorem [6l by considering the 
cases n odd and even separately. The only differences are that we replace 
r with p, the relation r^ + z"^ = 1 with z^ = p^ + 1, and the denominator 
(1 — kQj)^^'^ with (cqj- — 1)^/^, wherever it appears, where Cqj = —koj replaces 
koj. Unlike in S^, the case n even is satisfied for all admissible values of z. D 

Like in S^, the equilateral triangle presents particular interest, so let us say 
a bit more about it than in the general case of the regular n-gon. 

Corollary 3. Consider the 3-body with equal masses, m := rrii = m2 = m^, 
in H^. Then for any m > Q and z > 1, there are a positive and a negative uj 
that produce relative elliptic equilibria in which the bodies are at the vertices 
of an equilateral triangle that rotates in the plane z = constant. Moreover, for 
every u'^/m > there is a unique z > 1 as above. 

Proof. Substituting in system fl25]) a solution of the form 

(55) Xi = pcos{ujt + ai), yi = psm{ut + ai), Zi = z, 



with z = a/p^ + 1, «! = 0, a2 = 2tt/3, a^, = 47r/3, we are led to the equation 

8 LU^ 

56 -^ = — . 

^ ^ ^3(3^4 - 2^2 - 1)3/2 m 

The left hand side is positive for z > 1, tends to infinity when z ^ 1, and 
tends to zero when 2 — > oo. So for any z in this interval and m > 0, there are 
a positive and a negative u for which the above equation is satisfied. D 

As we already proved in the previous section, an equilateral triangle rotating 
in its own plane forms an elliptic relative equilibrium in S^ only if the three 
masses lying at its vertices are equal. The same result is true in H^, as we will 
further show. 

Proposition 3. In the 3-body problem in H^, if the bodies mi,m2,ms are 
initially at the vertices of an equilateral triangle in the plane z = constant for 
some z > 1, then there are initial velocities that lead to an elliptic relative 
equilibrium in which the triangle rotates in its own plane if and only if mi = 
7712 = m3. 
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Proof. The implication which shows that if mi = 7712 = rris, the rotating 
equilateral triangle is an elliptic relative equilibrium, follows from Theorem [TUl 
To prove the other implication, we substitute into equations (l28l) a solution of 
the form ([SS]) with i = 1, 2, 3, p := pi, P2, Ps, z := Zi = Z2 = Z3 = (p^ + 1)^/^, 
and «! = 0, Q!2 = 27r/3, «3 = in/S. The computations then lead to the system 

{mi + m,2 = C^'^ 
m2 + m3 = Cuj'^ 
m^ + mi = (uj^, 

where ( = v3(32^ — 2z'^ — 1)^/^/4. But for any z = constant with z > 1, the 
above system has a solution only for m-i = m,2 = m,^ = (uj'^/2. Therefore the 
masses must be equal. D 

The following result perfectly resembles Theorem [71 The proof works the 
same way, by just replacing the elliptical trigonometric functions with hyper- 
bolic ones and changing the signs to reflect the equations of motion in H^. 

Theorem 11. For all Lagrangian solutions in H^, the masses mi,m2 and 
m^ have to rotate on the same circle, whose plane must he orthogonal to the 
rotation axis, and therefore mi = 7712 = m,-^ . 

We will further prove an analogue of Theorem [HI 

Theorem 12. Consider the 3-body problem in H^ with equal masses, m : = 
mi = m2 = 777,3. Fix the body of mass mi at (0, 0, 1) and the bodies of masses 
m2 and m^ at the opposite ends of a diameter on the circle z = constant. 
Then, for any m > and z > 1, there are a positive and a negative uj, which 
produce elliptic relative equilibria that rotate around the z axis. 

Proof. Substituting into the equations of motion f[251) a solution of the form 

xi = 0, yi = 0, zi = 1, 

X2 = pcosci;t, 1/2 = p sin cut, Z2 = z, 

X3 = pcos(u;t + vr), 1/3 = psin(ci;t + vr), z^ = z, 

where p > and z > 1 are constants satisfying z"^ = p^ + 1, leads either to 
identities or to the algebraic equation 

.... 4^^ + l ^ ^ 

^ ' 4^3(^2 _ 1)3/2 ^• 

The function on the left hand side is positive for 2 > 1. Therefore, for every 
m > and z > 1, there are a positive and a negative u that lead to a geodesic 
elliptic relative equilibrium. The sign of u determines the sense of rotation. D 

Remark 6. For every u'^/m > 0, there is exactly one z > 1 that satisfies 
equation fISS]) (see Figure [S]). 
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Figure 5. The graph of the function f{z) 



423^2^1)3/2 for Z>1. 



6.3. Hyperbolic Relative Equilibria in H^. We now prove some results 
concerning hyperbolic relative equilibria. We first show that, in the n-body 
problem, hyperbolic relative equilibria do not exist along any given fixed ge- 
odesic of H^. In other words, the bodies cannot chase each other along a 
geodesic and maintain the same initial distances for all times. 

Theorem 13. Along any fixed geodesic, the n-body problem in H^ has no 
hyperbolic relative equilibria. 

Proof. Without loss of generality, we can prove this result for the geodesic 
X = 0. We will show that equations (1281) do not have solutions of the form 
(1^ with Xi = and (consequently) pi = 1, i = 1, . . . ,n. Substituting 



(59) 

into system 



Xi 



0, yi = sinh(ci;t 



Oii 



and Zi = cosh(ci;t + Oj) 



, the equation corresponding to the yi coordinate becomes 



(60) 






mj[smh.{ujt + aj) — cosh(aj 



a, 



sinh(u;t + ««)] 



sinh(aj — aj) 



a-i 



Assume now that «» > aj for all j ^ i. Let aM{i) be the maximum of all 
with j 7^ i. Then for t G {—aM{i)/^, —Oii/uj), we have that sinh(at + aj) < 
for all j 7^ i and sinh(at+aj) > 0. Therefore the left hand side of equation (l60l) 
is negative in this interval, so the identity cannot take place for all t G R. It 
follows that a necessary condition to satisfy equation (160|) is that aM{i) > ««• 
But this inequality must be verified for alii = 1, . . . ,n, a fact that can be 
written as: 

«! > a2 oi ai> a^ oi ... or ai > an, 

Ci2 > tti or a2> a-i or ... or 0^2 > an, 



O-n > Oil or an > ^2 or 



or an > a 



n-l- 
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The constants «!,...,«„ must satisfy one inequality from each of the above 
hues. But every possible choice implies the existence of at least one i and 
one j with i ^ j and ctj = aj. For those i and j, sinh(aj — aj) = 0, so 
equation (160|) is undefined, therefore equations (128|) cannot have solutions of 
the form (1591) . Consequently hyperbolic relative equilibria do not exist along 
the geodesic x = 0. D 

Theorem [13] raises the question whether hyperbolic relative equilibria do 
exist at all. For three equal masses, the answer is given by the following result, 
which shows that, in H^, three bodies can move along hyperbolas lying in 
parallel planes of R'^, maintaining the initial distances among themselves and 
remaining on the same geodesic (which rotates hyperbolically). The existence 
of such solutions is surprising. They rather resemble fighter planes flying in 
formation than celestial bodies moving under the action of gravity alone. 

Theorem 14. In the 3-body problem of equal masses, m := mi = m2 = m^, 
in H^, for any given m > and x ^ 0, there exist a positive and a negative u 
that lead to hyperbolic relative equilibria. 

Proof. We will show that qi(t) = {xi(t),yi{t), Zi{t)), i = 1,2, 3, is a hyperbolic 
relative equilibrium of system (l28l) for 

xi = 0, 2/1 = sinhut, z\ = cosh tut, 

X2 = X, y2 = psmh.u!t, Z2 = pcoshut, 

Xs = —X, y3 = psmh.ut, Zs = p cosh cut, 

where p = (1 + x^)^/^. Notice first that 

a;ia;2 + yiy2 - ziZ2 = xixs + yiy^ - ziz^ = -p, 

X2X2, + 2/22/3 - Z2Z-i = -2x^ - 1, 
•2|-2 -2 2 •2|-2 -2 •2|-2 -2 22 

Xi+Vi- Zi=uj , X2 + y2 - Z2 = x^ + y^ - z-^ = p UJ . 

Substituting the above coordinates and expressions into equations fl28l) . we are 
led either to identities or to the equation 

(a^ ) 4x^ + 5 ^ ^ 

^ ' 4x2|x|(a;2 + 1)3/2 ^' 

from which the statement of the theorem follows. D 

Remark 7. The left hand side of equation (1611) is undefined for a; = 0, but 
it tends to infinity when x —>■ and to when x — > ±oo. This means that 
for each uJ^ jrti > there are exactly one positive and one negative x (equal in 
absolute value), which satisfy the equation. 
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Remark 8. Theorem [H] is also true if, say, m := mi and M := 7712 = m^. 
Then the analogue of equation flMl) is 

m M p 



X^\X\(X' 



+ 1)1/2 4x2|x|(a;2 + 1)3/2 



and it is obvious that for any m, M > and x 7^ 0, there are a positive and 
negative cj satisfying the above equation. 

Remark 9. Theorem [H] also works for two bodies of equal masses, m : = 
nil = 1712, of coordinates 

Xi = —X2 = x,yi = y2 = p sinh ut, zi = Z2 = p cosh ut, 

where a; is a positive constant and p = (x^ + 1)^/^. Then the analogue of 
equation (16T]) is 

1 a;2 



) 



4x^1x1 (x^ + 1)3/2 j^ 
which obviously supports a statement similar to the one in Theorem [T 



6.4. Parabolic Relative Equilibria in H^. We now show that there are no 
parabolic relative equilibria. More precisely, we prove the following result. 

Theorem 15. The n-body problem in H^ has no parabolic relative equilibria. 

Proof. Let Xi,yi, and Zi be as in the definition of parabolic relative equilibria 
(^^. Then Xi = —hi + q, yi = ai + (q — bi)t, and a^ + (q — bi)t. Thus, the first 
component of the angular momentum is ^iTniai{bi — Ci) — ^iTni{bi — CiYt. 
It follows that bi = Ci because the first component of the angular momentum 
must be constant. But of + bf — c^ = —1, hence a^ = —1, which is impossible, 
since a,- is a real number. D 



7. Saari's conjecture 

In 1970, Don Saari conjectured that solutions of the classical n-body prob- 
lem with constant moment of inertia are relative equilibria, [SH], [OD]- The 
moment of inertia is defined in classical Newtonian celestial mechanics as 
1 J27=i ^iC[« ■ Qi! a function that gives a crude measure of the bodies' dis- 
tribution in space. But this definition makes little sense in S^ and H^ because 
Qi qi = =tl for every i = 1, . . . ,n. To avoid this problem, we adopt the 
standard point of view used in physics, and define the moment of inertia in S^ 
or H^ about the direction of the angular momentum. But while fixing an axis 
in S^ does not restrain generality, the symmetry of H^ makes us distinguish 
between two cases. 
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Indeed, in S^ we can assume that the rotation takes place around the z axis, 
and thus define the moment of inertia as 

n 

(62) \:=Y,mlx\^yl). 

1=1 

In H^, all possibilities can be reduced via suitable isometric transformations 
(see Appendix) to: (i) the symmetry about the z axis, when the moment of 
inertia takes the same form (l62l) . and (ii) the symmetry about the x axis, which 
corresponds to hyperbolic rotations, when — in agreement with the definition 
of the Lorentz product (see Appendix) — we define the moment of inertia as 

n 

(63) J:=5^m,(2/f-zf). 

The case of the parabolic roations will not be considered because there are no 
parabolic relative equilibria. 

These definitions allow us to formulate the following conjecture: 

Saari's Conjecture in S^ and H^. For the gravitational n-body problem in 
S^ and H^, every solution that has a constant moment of inertia about the 
direction of the angular momentum is either an elliptic relative equilibrium in 
S^ or H^, or a hyperbolic relative equilibrium in H^. 

By generalizing an idea we used in the Euclidean case, [23], [23], we can 
now settle this conjecture when the bodies undergo another constraint. More 
precisely, we will prove the following result. 

Theorem 16. For the gravitational n-body problem in S^ and H^, every so- 
lution with constant moment of inertia about the direction of the angular mo- 
mentum for which the bodies remain aligned along a geodesic that rotates el- 
liptically in S^ or H^, or hyperbolically in H^, is either an elliptic relative 
equilibrium in S^ or H^, or a hyperbolic relative equilibrium in H^. 

Proof. Let us first prove the case in which I is constant in S^ and H^, i.e. when 
the geodesic rotates elliptically. According to the above definition of I, we can 
assume without loss of generality that the geodesic passes through the point 
(0, 0, 1) and rotates about the 2;-axis with angular velocity uj(t) ^ 0. The 
angular momentum of each body is Lj = Tn^qj ® q^, so its derivative with 
respect to t takes the form 

with K = 1 in S^ and k = — 1 in H^. Since q^ Vq,t/K(q) = 0, it follows that 
Vq-f/K(q) is either zero or orthogonal to q^. (Recall that orthogonality here 
is meant in terms of the standard inner product because, both in S^ and H^, 
q^ Vq,f/,(q) = qi ■ Vq,t/,(q).) If Vq,t/,(q) = 0, then U = 0, so Lf = 0. 
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Assume now that Vq.f/K(q) is orthogonal to q^. Since all the particles are 
on a geodesic, their corresponding position vectors are in the same plane, 
therefore any linear combination of them is in this plane, so Vq.f/K(q) is in 
the same plane. Thus Vq-f/^lq) and q^ are in a plane orthogonal to the xy 
plane. It follows that Lj is parallel to the xy plane and orthogonal to the z 
axis. Thus the z component, LI, of Lj is 0, the same conclusion we obtained 
in the case Vq-f/K(q) = 0. Consequently, LI = Ci, where Cj is a constant. 

Let us also remark that since the angular momentum and angular velocity 
vectors are parallel to the z axis, L^ = IiUj{t), where I, = rrii^xj + yf) is the 
moment of inertia of the body mi about the 2;-axis. Since the total moment of 
inertia, I, is constant, and ui{t) is the same for all bodies because they belong 
to the same rotating geodesic, it follows that ^^^=1 '^i^{t) = Ii^(t) = c, where 
c is a constant. Consequently, cj is a constant vector. 

Moreover, since L^ = Ci, it follows that Iiu;(t) = q. Then every lj is 
constant, and so is every Zj, i = 1, . . . ,n. Hence each body of mass rrii has a 
constant 2;j-coordinate, and all bodies rotate with the same constant angular 
velocity around the 2;-axis, properties that agree with our definition of an 
elliptic relative equilibrium. 

We now prove the case J = constant, i.e. when the geodesic rotates hy- 
perbolically in H^. According to the definition of J, we can assume that the 
bodies are on a moving geodesic whose plane contains the x axis for all time 
and whose vertex slides along the geodesic hyperbola x = 0. (This moving ge- 
odesic hyperbola can be also visualized as the intersection between the sheet 
z > of the hyperboloid and the plane containing the x axis and rotating 
about it. For an instant, this plane also contains the z axis.) 

The angular momentum of each body is Lj = rrajqj Kl qj, so we can show 
as before that its derivative takes the form Lj = ?7ijqj Kl Vq.f/_i(q). Again, 
Vq,j?7_i(q) is either zero or orthogonal to qj. In the former case we can draw 
the same conclusion as earlier, that Lj = 0, so in particular L^ = 0. In the 
latter case, qj and Vq-t/_i(q) are in the plane of the moving hyperbola, so 
their cross product, qj Kl Vq.f/_i(q) (which differs from the standard cross 
product only by its opposite z component), is orthogonal to the x axis, and 
therefore L^ = 0. Thus L^ = in either case. 

From here the proof proceeds as before by replacing I with J and the z axis 
with the X axis, and noticing that L^ = Jju;(t), to show that every rrii has a 
constant Xi coordinate. In other words, each body is moving along a (in general 
non-geodesic) hyperbola given by the intersection of the hyperboloid with a 
plane orthogonal to the x axis. These facts in combination with the sliding of 
the moving geodesic hyperbola along the fixed geodesic hyperbola x = are 
in agreement with our definition of a hyperbolic relative equilibrium. D 
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8. Appendix 

8.1. The Weierstrass model. Since the Weierstrass model of the hyperbolic 
(Bolyai-Lobachevsky) plane is little known, we will present here its basic prop- 
erties. This model appeals for at least two reasons: (i) it allows an obvious 
comparison with the sphere, both from the geometric and analytic point of 
view; (ii) it emphasizes the differences between the Bolyai-Lobachevsky and 
the Euclidean plane as clearly as the well-known differences between the Eu- 
clidean plane and the sphere. As far as we are concerned, this model was the 
key for obtaining the results we proved for the n-body problem for n < 0. 
The Weierstrass model is constructed on one of the sheets of the hyperboloid 



x'^ + y'^ — z"^ = —1 in the 3-dimensional Minkowski space Ai^ := (R'^, H), in 



which a □ b = Uxbx + CLyby — Uzbz, with a 



yClxy Cly, 



ttz) and b = (6^, 



by,bz 



represents the Lorentz inner product. We choose the 2; > sheet of the 
hyperboloid, which we identify with the Bolyai-Lobachevsky plane H^. 

A linear transformation T: A4'^ — *> Ai^ is orthogonal if T(a) E]T(a) = aEla 
for any a G Ai^. The set of these transformations, together with the Lorentz 
inner product, forms the orthogonal group 0{Ai^), given by matrices of de- 
terminant ±1. Therefore the group SO{At^) of orthogonal transformations 
of determinant 1 is a subgroup of 0{Ai^). Another subgroup of 0{Ai^) 
is G{At^), which is formed by the transformations T that leave H^ invari- 
ant. Furthermore, G{Ai^) has the closed Lorentz subgroup, Lor(A^^) : = 
G{M^)nSO{M'^). 

An important result is the Principal Axis Theorem for Lor(AI^), [26], [M] . 
Let us define the Lorentzian rotations about an axis as the 1-parameter sub- 
groups of Lor(A^'^) that leave the axis pointwise fixed. Then the Principal Axis 
Theorem states that every Lorentzian transformation has one of the forms: 



A = P 



A 



cos 6 

sin^ 





1 
cosh s 
sinh s 



- sin 6 


0" 


cos 6 








1 




sinhs 
coshs 



P- 



P- 



or 



1 -t t 

A = P t 1 - tV2 tV2 P- 

t -tV2 1 + 1^/2 

where 6 G [0,27r), s,t G R, and P G LoT{Ai^). These transformations are 
called elliptic, hyperbolic, and parabolic, respectively. The elliptic transfor- 
mations are rotations about a timelike axis — the z axis in our case; hyperbolic 
rotations are rotations about a spacelike axis — the x axis in our context; and 
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parabolic transformations are rotations about a lightlike (or null) axis, repre- 
sented here by the line x = 0, y = z. This result resembles Euler's Principal 
Axis Theorem, which states that any element of 5*0(3) can be written, in some 
orthonormal basis, as a rotation about the z axis. 

The geodesies of H^ are the hyperbolas obtained by intersecting the hy- 
perboloid with planes passing through the origin of the coordinate system. 
For any two distinct points a and b of H^, there is a unique geodesic that 
connects them, and the distance between these points is given by d{sL, b) = 
cosh"^(-aHb). 

In the framework of Weierstrass's model, the parallels' postulate of hyper- 
bolic geometry can be translated as follows. Take a geodesic 7, i.e. a hyperbola 
obtained by intersecting a plane through the origin, O, of the coordinate sys- 
tem with the upper sheet, 2; > 0, of the hyperboloid. This hyperbola has two 
asymptotes in its plane: the straight lines a and b, intersecting at O. Take a 
point, P, on the upper sheet of the hyperboloid but not on the chosen hyper- 
bola. The plane aP produces the geodesic hyperbola a, whereas bP produces 
p. These two hyperbolas intersect at P. Then a and 7 are parallel geodesies 
meeting at infinity along a, while /5 and 7 are parallel geodesies meeting at in- 
finity along b. All the hyperbolas between a and (3 (also obtained from planes 
through O) are non-secant with 7. 

Like the Euclidean plane, the abstract Bolyai-Lobachevsky plane has no 
privileged points or geodesies. But the Weierstrass model has some convenient 
points and geodesies, such as the point (0, 0, 1) and the geodesies passing 
through it. The elements of Lor(AI^) allow us to move the geodesies of H^ 
to convenient positions, a property we frequently use in this paper to simplify 
our arguments. Other properties of the Weierstrass model can be found in [29] 
and [58]. The Lorentz group is treated in some detail in [2], but the Principal 
Axis Theorems for the Lorentz group contained in [2] and [58] fails to include 
parabolic rotations, and is therefore incomplete. 



8.2. History of the model. The first researcher who mentioned Karl Weier- 
strass in connection with the hyperboloidal model of the Bolyai-Lobachevsky 
plane was Wilhelm Killing. In a paper published in 1880, [ID], he used what 
he called Weierstrass's coordinates to describe the "exterior hyperbolic plane" 
as an "ideal region" of the Bolyai-Lobachevsky plane. In 1885, he added 
that Weierstrass had introduced these coordinates, in combination with "nu- 
merous applications," during a seminar held in 1872, [12], pp. 258-259. We 
found no evidence of any written account of the hyperboloidal model for the 
Bolyai-Lobachevsky plane prior to the one Killing gave in a paragraph of [12] , 
p. 260. His remarks might have inspired Richard Faber to name this model 
after Weierstrass and to dedicate a chapter to it in [29], pp. 247-278. 
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